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Abstract. 

Ice crystals in clouds exist in a virtually limitless variation of geometries. The most basic shapes of ice crystals are columnar 
or plate-like hexagonal prisms with aspect ratios determined by relative humidity and temperature. However, crystals in ice 
clouds generally display more complex structures owing to aggregation, riming and growth histories through varying tempera- 
ture and humidity regimes. Crystal shape is relevant for cloud evolution as it affects microphysical properties such as fall speeds 
and aggregation efficiency. Furthermore, the scattering properties of ice crystals are affected by their general shape, as well 
as by microscopic features such as surface roughness, impurities and internal structure. To improve the representation of ice 
clouds in climate models, increased understanding of the global variation of crystal shape and how it relates to, e.g., location, 
cloud temperature and atmospheric state is crucial. Here, the remote sensing of ice crystal macroscale and microscale structure 
from airborne and space-based lidar depolarization observations and multi-directional measurements of total and polarized 
reflectances is reviewed. In addition, a brief overview is given of in situ and laboratory observations of ice crystal shape as well 
as the optical properties of ice crystals that serve as foundations for the remote sensing approaches. Lidar depolarization is gen- 
erally found to increase with increasing cloud height and to vary with latitude. Although this variation is generally linked to the 
variation of ice crystal shape, the interpretation of the depolarization remains largely qualitative and more research is needed 
before quantitative conclusions about ice shape can be deduced. The angular variation of total and polarized reflectances of ice 
clouds has been analyzed by numerous studies in order to infer information about ice crystal shapes from them. From these 
studies it is apparent that pristine crystals with smooth surfaces are generally inconsistent with the data and thus crystal impu- 
rity, distortion or surface roughness is prevalent. However, conclusions about the dominating ice shapes are often inconclusive 
and contradictory and are highly dependent on the limited selection of shapes included in the investigations. Since ice crystal 
optical properties are mostly determined by the aspect ratios of the crystal components and their microscale structure, it is 
advised that remote sensing applications focus on the variation of these ice shape characteristics, rather than on the macroscale 
shape or habit. Recent studies use databases with nearly continuous ranges of crystal component aspect ratio and/or roughness 
levels to infer the variation of ice crystal shape from satellite and airborne remote sensing measurements. Here, the rationale 
and results of varying strategies for the remote sensing of ice crystal shape are reviewed. Observed systematic variations of ice 
crystal geometry with location, cloud height and atmospheric state suggested by the data are discussed. Finally, a prospective 


is given on the future of the remote sensing of ice cloud particle shapes. 
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1 Introduction 


Archaeological findings reveal a deeply rooted human obsession for atmospheric optical phenomena such as rainbows and 
ice cloud halos. Our pre-historic ancestors likely attached divine or meteorological meanings to such optical phenomena as 
suggested by ancient rock drawings (Ping-Yii and Needham, 1959; Sassen, 1991, 1994). Detailed observations of rainbows 
and ice cloud halos were later described by Aristotle and other scholars of antiquity. It took until the age of enlightenment, 
however, for the first qualitatively correct interpretations of these phenomena to be documented (Nussenzveig, 1977; Greenler, 
1990; Tape, 1994; Tape and Moilanen, 2006). For example, it was the French physicist and priest Edme Marriotte who, in 
1691, was probably the first to correctly attribute the commonly observed halo at an angle of 22° around the sun to randomly 
oriented hexagonal ice crystals, although his ideas did not gain much attention until being revived 116 years later by Thomas 
Young (Tape and Moilanen, 2006). It was later realized that this hexagonal ice model also offers explanations of parhelia 
(Csundogs’) and other arcs, when crystals are assumed to be oriented. One can thus argue that Mariotte’s interpretation of the 
halo observations was the first successful remote sensing of ice crystal shape. 

While liquid drops generally take shape of virtually perfect spheres, the observed variety of cloud ice crystal shape is seem- 
ingly limitless (Baran, 2009). Under most atmospheric conditions and undisturbed, water ice grows as hexagonal prisms owing 
to its fundamental hexagonal molecular structure (generally denoted as J;,), where the growth rate of the basal and prismatic 
planes of the prisms are determined by temperature and humidity (Pauling, 1935; Bailey et al., 2012; Harrington et al., 2013). 
However, because of varying atmospheric conditions throughout crystal evolution, aggregation, fracturing, riming, and other 
stochastic effects, observed cloud ice crystals are generally far more complex than simple hexagonal prisms (e.g., Baran, 2009; 
Bailey and Hallett, 2009; Bailey et al., 2012). Moreover, recent advances in electron microscope technology allows imaging 
of the microstructure of ice surfaces, which reveals that roughness structures are generally prevalent (Pfalzgraff et al., 2010; 
Neshyba et al., 2013; Magee et al., 2014). Mixing of the hexagonal ice molecular structure with the cubical structure (de- 
noted as J.) can lead to such microscale roughness structures on the surfaces, in addition to more macroscale deviation 
of ice crystals from their fundamental hexagonal form (Kuhs et al., 2012; Malkin et al., 2012, 2015; Murray et al., 2015; 
Hudait and Molinero, 2016). Furthermore, sublimation and riming changes the surfaces of ice crystals (Ono, 1969; Avila et al., 


2009; Pfalzgraff et al., 2010; Magee et al., 2014). 


The shape and size of cloud particles greatly affects their radiative and microphysical properties (Yang et al., 2015; Furtado et al., 


2015; Fridlind et al., 2016; Baran et al., 2016; Liou and Yang, 2016). The omnipresence of clouds on Earth means that seem- 


ingly small changes in cloud particle properties can have profound effects on the radiation balance and precipitation (Stephens et al., 


1990). In turn, cloud radiative properties affect the evolution of clouds owing to interactions between radiation and cloud micro- 
physics (e.g., Gu and Liou, 2000; Gu et al., 2011; Russotto et al., 2016; Baran et al., 2016). While the dependences of radiative 
and microphysical properties on particle size are fairly well understood, knowledge on how they are affected by the shape 
of ice crystals is still limited. The current incomplete knowledge on the natural variation of ice crystal shape contributes to 
our inability to reliably represent ice clouds in global circulation models (GCMs). Biases in ice optical properties can lead 


to direct biases in modeled ice cloud radiative effects (Yi et al., 2013; van Diedenhoven et al., 2014a; Liou and Yang, 2016). 
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Furthermore, the relatively large uncertainties of cloud physics and the large influence of clouds on the radiation balance 
means that weakly constrained variables in GCM cloud physics parameterizations are generally adjusted to bring the simu- 
lated global long- and shortwave radiation in agreement with measurements (e.g., Schmidt et al., 2014). Through this process, 
any biases in ice particle optical properties would thus be falsely compensated by inadvertently introducing biases in cloud 
physical parameters. Moreover, biases in cloud particle shape and related microphysical quantities such as fall speeds can lead 
to biases in modeled cloud properties such as thermodynamic phase, precipitation and cloud fraction (Furtado et al., 2015). 
A new paradigm in ice microphysical modeling is emerging in which ice crystal shape characteristics are explicitly predicted 
(Hashino et al., 2007, 2011; Harrington et al., 2013), although a lack of observations for evaluation leave these models largely 
unconstrained. In addition to ice crystal shape, global constrains on the ice crystal size are also crucial in order to evaluate 
and improve the representation of clouds in models. However, the large uncertainty on the variation of ice particle shape, and 
thus optical properties, in turn lead to substantial uncertainties in ice particle size retrieved using satellite-based shortwave 
infrared measurements, such as by the Moderate Resolution Imaging Spectroradiometer (MODIS) (van Diedenhoven et al., 
2014b; Holz et al., 2016). In situ and laboratory measurements of the past decades have lead to an increasingly detailed picture 
of the macro- and microscale structure of ice under varying conditions and nucleation processes (Bailey and Hallett, 2009). 
However, this picture is incomplete and likely biased, given the many technical and practical limitations of in situ and labora- 
tory measurements and the uncertain and complex variation of ice nucleation and growth conditions and processes occurring in 
Earth’s ice-containing clouds (cf. Bailey et al., 2012). Remote sensing observations of ice shape on global and regional scales 
are crucial for the evaluation and refinement of parameterizations of ice crystal nucleation and evolution derived from in situ 
and laboratory observations and theory. 

This chapter offers a review of the current state of remote sensing of ice crystal macro- and microscale structure. We focus 
on remote sensing techniques based on passive total and polarized reflectance measurements. In addition, the potential of lidar 
with polarization capability to infer information about ice crystal shape is reviewed. A discussion on polarimetric radar is 
excluded here since current polarimetric radars are mostly sensitive to precipitation rather than cloud ice, and much of the 
polarimetric radar signal is related to particle orientation rather than shape (Miao et al., 2003). First, in section 2 we will briefly 
summarize the in situ and laboratory observations that serve as foundations of the retrieval approaches. Since by definition 
remote sensing necessitates the interpretation of measured electromagnetic radiation, we will review the optical properties of 
ice crystals in section 3. Active and passive techniques for the remote sensing of ice crystals shape are reviewed in section 4. 


A prospective will be given in section 5, before concluding the chapter in section 6. 


2 Insitu and laboratory measurements 


In 1665, English natural philosopher Robert Hooke included some of the first detailed drawings of snow crystals in his book 
Micrographia (Hooke, 1665). The first detailed photographs were collected around the turn of the 20th century by enthusiasts 
such as Wilson Bentley who documented many snow crystals of various common and uncommon shapes (Bentley, 1927). The 


first known in situ aircraft observations of ice crystal shape were performed during the second world war by German meteorol- 
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ogist Ludwig Weickmann, who performed about 100 meteorological flights in an open cockpit airplane holding lacquer-coated 
shingles to capture and replicate ice crystals. Weickmann published results and their interpretation after the end of the war 
(Weickmann, 1945; aufm Kampe et al., 1951). To date, many in situ and laboratory observations of ice crystal shape have 
been conducted. Here, we will give a brief overview of the general and detailed shape characteristics measured in situ and in 


laboratories. Such observations are essential to direct and constrain remote sensing applications and to interpret their results. 
2.1 General classification of ice habits 


Crystal shapes are commonly classified according to the scheme of Magono and Lee (1966) that contains 80 classes of snow 
and ice crystals. This scheme in turn is based on previous work by Nakaya (1954), while a proposed expansion of the clas- 
sification was recently offered by Kikuchi et al. (2013). The identified ice crystal shapes are commonly referred to as ice 
*habits’. Nakaya (1954) and Magono and Lee (1966) also showed that the occurrence of certain habits could be linked to the 
meteorological conditions in which they grew and they provided diagrams of occurrence of certain ice habits as a function 
of temperature and humidity. Around the same time, laboratory experiments in cloud chambers were performed that success- 
fully reproduced many of the observed transitions of habits at certain temperatures and humidities (aufm Kampe et al., 1951; 
Mason, 1953; Hallett and Mason, 1958; Mason et al., 1963). More recently, the most comprehensive habit diagram for atmo- 
spheric ice crystals was published by Bailey and Hallett (2009), which is shown in Fig. 1. As is apparent from Fig. 1, there 
are two major differentiating characteristics of ice crystals: (1) plate-like versus columnar structures and (2) single crystals 
versus polycrystalline particles. Plate-like crystals form when the basal planes grow faster than the prismatic facets, while the 
opposite relation of facet growth rates lead to columnar crystals. Polycrystals can form for various reasons, for example, when 
several crystal components grow relatively independently from a common core, when crystals aggregate, or when conditions 
vary during crystal growth. Comparing habits in different cloud types, one general observation is that compact and aggre- 
gated ice crystals with plate-like components often occur in tropical deep convection (e.g. Noel et al., 2004; Connolly et al., 
2005; Um and McFarquhar, 2009; Baran, 2009), while bullet rosettes are often found in non-convective ("in situ") cirrus (e.g., 
Lawson et al., 2010; Fridlind et al., 2016). 

The diagram produced by Bailey and Hallett (2009) (Fig. 1) is partly based on images of individual ice crystals taken in 
clouds by the Cloud Particle Imager (CPI) aircraft probe (Lawson et al., 2001, 2006). The highly successful CPI probe was 
designed in 1997 by SPEC inc. and can image and count particles in the size range of 15-2500 wm, with a nominal resolution 
of 2.3—jum. As discussed by Bailey and Hallett (2009), CPI images reveal many complex, polycrystalline shapes that were not 
represented by laboratory studies before. Automatic habit classification applied to the images generally yield a dominance 
of unclassified or irregular’ crystals. Recently developed, more advanced classification schemes (e.g., Lindqvist et al., 2012) 
add more classes to previous CPI classification programs, but still yield an abundance of unclassified crystals. The general 
dominance of crystals classified as ’irregular’ may be partly attributed to the issues with the automatic classification of 2D 
images (Stoelinga et al., 2007). Moreover, the CPI resolution is generally insufficient to classify small crystals into classes other 
than "small irregular" or ’quasi-spherical’. Although true quasi-spherical crystals have been observed (Jarvinen et al., 2016), 


for most cases the term ’quasi-spherical’ is arguably a misnomer, since on closer inspection these crystals generally appear 


0.6 +—elamnar o'verystalline Regime itelike > Single Crystals 


*+¢—- approximate maximum 


hollow ended in situ ice supersaturation 


& columns 


0.5 -long bullets and columns ‘-. 1 . due to fall and ventilation 
5 needle assemblages 
= rosettes bullet rosettes with > Of large plates. 
% long solid bullets Send side planes - 
z * a 
long columns crossed E a 
2 plates © : 
rT) d plates 6 3: 
+4 needle solid bullet rosettes crossed p' t o- 
3 rosettes side planes © = 6 
” solid columns spearheads . e& E 
0 gohei twins 3 a 5 
2 column/plate scrolls tie r-) 
combinations combinations assemblages thin wo 
of irregulars of thin plates skeletal $ & 
plates “ =. 
° 
£ 


short solid columns compact assemblages thin plates 


of plates, thick & thin plates thick plates 
columns 


———————Polycrystalline Re gime———__—__> Single Crystals 


columnar plate-like 


Ice Supersaturation 


eve eg dh ¥ FH 


-70 -60 -50 -40 -30 -20 -10 0 
Temperature (°C) 


Figure 1. Diagram showing the general occurrence of ice crystal habits at varying temperatures and water vapor supersaturations 
with respect to ice. Habit names are in the top diagram, while example images taken by the CPI instrument are in the lower 


diagram. See Bailey and Hallett (2009) for more details. Figure reproduced from Bailey and Hallett (2009). 
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as compact polyhedra with many small facets (Gonda and Yamazaki, 1978; Bailey and Hallett, 2009) or irregular shapes with 
bulges and protuberances (Nousiainen et al., 2011; Lopez and Avila, 2012). 

In summary, in situ observations are found to follow the habit diagrams as a function of temperature and humidity, but 
seemingly only for crystals that have grown under relatively static conditions. Given the variability of atmospheric humidity, 
the fact that growing ice tends to fall to warmer temperatures and the presence of up- and downdrafts in clouds, static growth 
conditions are an exception rather than the rule, especially in convectively-driven clouds. Furthermore, processes such as crystal 


aggregation, fracturing and riming add more complexity to natural crystals. 
2.2 Aspect ratios of ice crystals 


Traditionally, ice crystal classifications and observations have focused on the large-scale appearance of the crystals, that is, on 
the habits. As will be discussed in section 3, smaller scale details of the crystals are also of importance for their optical and 
microphysical properties, if not more important than habit. One such property is the relation between the thickness and length 
of particles or components of polycrystals. For a hexagonal prism, the fraction between crystal basal plane width W and prism 


length L is commonly referred to as aspect ratio a, which can be defined as 
a= L/W. (1) 


In this definition the aspect ratio of a column is greater than unity and that of a plate smaller than unity. However, this choice is 
arbitrary and confusingly some studies define aspect ratio as the inverse of Eq. 1. Alternatively, van Diedenhoven et al. (2016a) 


proposed a definition that limits the aspect ratio to below unity for both plates and columns, viz. 

ac. = erat (2) 
When this definition is used, it needs to be specified separately whether crystals are column- or plate-like. Also note that 
the term ’aspect ratio’ is sometimes used for the ratio between two perpendicular dimensions of complex ice crystals (e.g., 
Korolev and Isaac, 2003) instead of the ratio between dimensions of components of those complex crystals, as the definition 
used here. 

Several early studies on dimensional relationships of ice crystals were reported on by, e.g., Mason (1953), Ono (1969) and 
Auer and Veal (1970) (See Um et al. (2015) for a completer list). Auer and Veal (1970) deduced power-law relationships be- 
tween the orthogonal dimensions of relatively simple particles such as plates, columns and dendrites that are still much used 
today as references for optical property calculations (e.g., Yang et al., 2015). Dimensions of more complex particles, such as 
bullet rosettes and aggregates of plates, were estimated by, e.g., Heymsfield and Andrew (1972), Mitchell and Arnott (1994), 
Xie et al. (2011), Um and McFarquhar (2007, 2009); Um et al. (2015) and Fridlind et al. (2016). However, most of these es- 
timates are based on very limited data and unconstrained assumptions and different techniques can lead to conflicting results 
(Um et al., 2015). For plates, reported aspect ratios range from 0.01 to 0.5, while aspect ratios of columns generally range 
from 1.5 to 5 (Auer and Veal, 1970; Um et al., 2015). In general, aspect ratios increasingly deviate from unity as the maximum 


dimension increases. For bullet rosettes and their aggregates, reported aspect ratios of the arms also range from 1.5 to 5.0 and 
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are found to decrease towards unity with decreasing number of bullets attached to the rosettes (Um et al., 2015; Fridlind et al., 
2016). Um et al. (2015) found that the aspect ratios of plates, columns and bullets in rosettes are not or only weakly sen- 
sitive to temperature, although a substantial temperature dependence may be expected from theory (Chen and Lamb, 1994; 
Harrington et al., 2013). Note that the crystals for which the dimensions can be determined generally comprise only a very 
limited subset of crystals in the observed cloud. Furthermore, the values listed above are mostly derived by the functional 
relationships between crystal (component) length and width that are determined from the data by the individual studies. These 
relationships, however, do not reflect information about the spread of data from which they are derived, which can be consid- 
erable (e.g., Um et al., 2015; Fridlind et al., 2016). For other complex crystal habits, such as aggregates of plates or columns, 
determining the dimensional relationships of their components from imagery is more problematic causing their aspect ratios to 


be more uncertain (e.g., Um and McFarquhar, 2009). 
2.3 Microscale structure of ice crystals 


A perfect, pristine’, ice crystal exclusively consisting of hexagonal molecular structure (J;,) will grow smooth basal and 
prismatic surfaces. However, it has been long known that imperfections in the J; mesh of growing water ice are prevalent, 
even under stable growth conditions (Pauling, 1935; Mason et al., 1963; Cross, 1969; Mizuno and Yukiko, 1978). Electron 
microscopic observations of ice crystals grown in the laboratory under different conditions reveal varying imperfections or 
roughness structures on the ice surfaces (Pfalzgraff et al., 2010; Sazaki et al., 2010; Neshyba et al., 2013; Magee et al., 2014). 
As discussed in section 3.2, such surface distortions strongly affect the scattering properties of ice crystals and need to be con- 
sidered in techniques for the remote sensing of ice crystal shapes. Ice crystal surface roughness generally decreases asymmetry 
parameters and depresses halos in phase functions at 22° and 46° (Yang et al., 2008a; Neshyba et al., 2013; van Diedenhoven, 
2014; van Diedenhoven et al., 2014a; Yang et al., 2015). 

An example of roughness structures observed on laboratory-grown ice crystals under various conditions is shown in Fig. 2 
(Magee et al., 2014). Magee et al. (2014) found that such roughness structures are prevalent and not confined to a narrow range 
of macroscopic morphology, substrate, temperature, humidity, or growth rates under most conditions. Although they did not 
detect a systematic dependence on the degree of supersaturation or the rate of growth, clear surface morphology differences 
between growing and sublimating crystals were observed. Specifically, sublimating crystals develop scalloped depressions and 
sharp ridges, while more linear structures were observed for growing crystals (see Fig. 2, also cf. Pfalzgraff et al., 2010). 
Indirect observations of roughness are provided by the Small Ice Detector-3 (SID-3) probe mounted in cloud chambers and 
on aircraft. This probe measures the light scattered in near-forward direction by single ice crystals. As discussed also in 
section 3.2, roughness structures will randomize the scattering angles at which light is scattered forward resulting in ’speckled’ 
images measurement by the SID-3 probe, while pristine particles lead to more organized, symmetric measurements patterns. A 
subjective, arguably arbitrary roughness parameter is deduced from these measurements to quantify the degree of roughness. 
Using the SID-3 probe mounted on aircraft sampling cirrus and convective clouds, Ulanowski et al. (2014) and Schmitt et al. 
(2016), respectively, found that roughness is ubiquitous. Somewhat in contrast with the results of Neshyba et al. (2013) and 


Magee et al. (2014), cloud chamber experiments with the SID-3 measurements conducted by Schnaiter et al. (2016) indicate a 
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Figure 2. High-magnification images of hexagonal ice crystals acquired by environmental scanning electron microscopy reveal- 
ing roughness structures on the crystal surfaces. The left image shows the crystal under vapor growth conditions. The center 
and right images show the crystal at the initial and a progressed stage of sublimation, respectively. For details see Magee et al. 


(2014). Images courtesy of Dr. Nathan Magee at The College of New Jersey. 


clear correlation between the small-scale crystal complexity and the volume mixing ratio of available condensable water vapor 
determining the growth rate. Furthermore, Ulanowski et al. (2014) reported on similar roughness in the growth and sublimation 
zones of cirrus observed by the SID-3 probe, which also seems to contradict conclusions from electron microscope observations 
(Magee et al., 2014). Clearly, more research needs to be conducted to investigate the presence and level of roughness in the 
varying types of ice clouds and under varying conditions. 

Although roughness structures are found to be prevalent on natural ice crystals, this appears to be in conflict with the frequent 
sighting of halos around the sun at 22° that are attributable to hexagonal crystals with smooth surfaces (Sassen et al., 2003; 
Tape and Moilanen, 2006; Verschure, 1998). For example, Sassen et al. (2003) reported that 37.3% of the daytime ice cloud 
observations in their ~10 year record over Salt Lake City, Utah, showed indications of the 22° halo, with bright and prolonged 
halos occurring in 6% of the record. van Diedenhoven (2014) reconciled the abundance of rough ice surfaces detected with the 
prevelance of halos by considering cirrus clouds to contain mixtures of crystals with varying roughness levels. van Diedenhoven 
(2014) concluded that that the contribution by pristine crystals to the total scattering cross section needs to be greater than only 
about 10% in the case of compact particles or columns, and greater than about 40% for plates for the 22° halo feature to 
be present in scattering phase functions. These results indicate that frequent sightings of 22° halos are not inconsistent with 
the observed dominance of rough ice crystals. From SID-3 measurements in a halo producing cloud, Schmitt et al. (2016) 
also concluded that it does not require high concentrations of halo producing particles to produce an observable halo. The 
occurrence of 22° halo and other optical phenomena attributed to pristine crystals in some clouds but the lack there-of in other 
clouds illustrates that ice crystal microscale morphology is highly variable. Note that the contrast and brightness of halos also 


depends on the cloud optical thickness (Kokhanovsky, 2008). 
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Roughness structures can be effectively observed by present-day instruments, but the quantification of roughness remains 
elusive. One metric to quantify roughness from electron micoscope imagery was proposed by Neshyba et al. (2013) and rela- 
tionships between this roughness metric and roughness parameterizations used for calculations of optical properties of rough 
ice crystals, as discussed in section 3.2, were also provided. Unfortunately, Magee et al. (2014) found that this metric proposed 
by Neshyba et al. (2013) is dependent on image magnification and resolution and is thus rather subjective. As surface rough- 
ness has profound effects on the scattering properties of ice crystals, and an increasing number of studies of the formation 
roughness structures under varying conditions are conducted, it is advisable that a more universal definition of roughness ought 
to be developed. 

Apart from roughness structures on crystal surfaces, other small-scale ice crystal imperfections occur. Specifically, hollow 
endings in ice columns or bullets are frequently reported (Heymsfield et al., 2002; Schmitt et al., 2007; Bailey and Hallett, 
2009; Smith et al., 2015). Bailey et al. (2004) found that hollowness appears at relatively high supersaturations of 50% and the 
depth of hollowness increases with increasing supersaturation. The hollowness generally appears to be conical (Heymsfield et al. 
2002), but step-wise hollowness is also reported (Smith et al., 2015). Other observed crystals imperfections include air bubbles 


and soot-inclusions (e.g., Hong and Minnis, 2015; Panetta et al., 2016). 


3 Ice crystal optical properties 


For the interpretation of ice cloud remote sensing observations, appropriate optical properties of ice crystals are essential. Be- 
cause ice crystals have non-spherical, complex structures and are generally much larger than solar wavelengths, exact methods 
for the calculation of optical properties, such as Lorenz—Mie theory (Mie, 1908) or the T-matrix approach (Mishchenko, 1991), 
are not suitable. Commonly, ice optical properties are approximated using geometric optics (GO) calculations, sometimes in 
combination with other methods. Already in the 17th century, GO principles were applied by Edme Mariotte and others to 
explain ice cloud halos. Modern applications are based on Monte Carlo ray tracing techniques that combine calculated paths 
of a large, finite number of light rays through the particle of interest assuming random orientations of the incoming ray. Early 
implementations of ray tracing were presented by Jacobowitz (1971) and Wendling et al. (1979), but more versatile computer 
codes that can be applied to complex particles and take into account light polarization are provided by Takano and Liou (1989) 
and Macke (1993) and are still much used today with some modifications. Although conventional GO has several shortcomings, 
mostly these issues are limited to its application to small and pristine particles with smooth surfaces (Bi et al., 2014; Yang et al., 
2015). A comprehensive review of optical property calculations using conventional and modified geometric optics and other 
methods is given by Yang et al. (2015). Here, we aim to summarize the optical properties that are needed for the simulation of 


remote sensing data and to give a brief review of the dependency of optical properties on particle shape characteristics. 
3.1 Definitions of optical properties 


In order to solve the radiative transfer problem in the case of an ice cloud, the cloud layer optical thickness, single scattering 


albedo and phase matrix are needed (van de Hulst, 1957). The dimensionless optical thickness 7 of a uniform ice cloud layer 
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of thickness AZ is given by 
tT =AZ Nrot Te (3) 


where Ni.¢ 1s the total number concentration of ice crystals and oe is the mean ice crystal extinction cross section. The 
extinction cross section for individual crystals can be written as a combination of absorption a, and scattering a, cross sections, 


Le. 
Ce =O04¢ +Os: (4) 


The fraction of light that is scattered rather than absorbed is given by the single scattering albedo w defined as 


Ba (5) 
Oe 

The extinction cross section of an ice crystal is largely determined by the projected area A, of the ice crystal, i-e., 

Oe = Qe Ap, (6) 


where Q, is the extinction efficiency. The projected area is defined as the crystal’s average projection when an infinite number 
of random orientations are applied (Macke, 1993). For convex particles, A, is equal to the total surface area divided by 4 (Vouk, 
1948), while for particles with concave parts the projected area can be calculated using a Monte Carlo approach (Macke, 1993). 
In the conventional GO approximation, @. is 2 by definition (Berg et al., 2011), but actual values depend on particle size and 
shape (Yang et al., 2015). For particles with sizes comparable to the considered wavelength, ray interference leads to oscilla- 
tions of Q, around 2 that decrease in strength with increasing particle size. For particles much smaller than the wavelength, Q, 
strongly decreases as it enters the Rayleigh scattering regime. Note however that generally the extinction coefficient is not very 
relevant for most remote sensing applications as the optical thickness is generally retrieved directly. However, the extinction 
coefficient is important for radiative transfer calculations based on model results (e.g., van Diedenhoven et al., 2012b; Bi et al., 
2014; Baran et al., 2016). 

The intensity and polarization state of light can be described by the Stokes vector S = (1,Q,U,V), where I is the total 
intensity, Q and U describe the linear polarization state and V describes the circular polarization state (van de Hulst, 1957). 
For each scattering event on randomly oriented ice crystals, the alteration of propagation direction and polarization state of 


incoming light is described by the symmetric phase matrix P, which can be defined through (van de Hulst, 1957). 
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where subscript "in" denotes the incoming Stokes vector, while "s" denotes the Stokes vector after scattering. Furthermore, r 


is the distance between the particle and the location at which Sg is evaluated. Note that all elements are a function of scattering 
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angle O. The first element of the phase matrix Pj; is referred to as the scattering phase function and is the only element to 
be considered when polarization is ignored. The other elements describe the alteration of the light’s polarization state at each 
scattering event. For example, the Pz element determines the linear polarization of light after one scattering event for incoming 
unpolarized light. The degree of linear polarization (DoLP) is given by P,2/P,;. Circular polarization can usually be ignored 
(Stamnes et al., 2016), but is relevant for lidars emitting circular polarized light (e.g., van Diedenhoven et al., 2009). 

In the geometrics optics regime, the total scattering phase function includes contributions from refraction into and out of the 


particle and internal reflections and contributions from diffraction (Macke et al., 1996b), viz. 
1 
P11(0) = om [(2u —1)Prr(9) + Pair(9)], (8) 


where Pai¢ and Pry are the contributions by diffraction and internal refraction and reflections, respectively. For radiative 
transfer applications that aim to approximate the upward and downward fluxes of light, e.g., so-called two-stream applications 
(Coakley and Chylek, 1975), the full phase function is not needed and only its first moment, the asymmetry parameter g, is 


considered. The asymmetry parameter is defined as 
g= J Pol) cos(@) sin(@)dO. (9) 
0 


The reflection of an ice cloud layer at a given wavelength and scattering geometry is primarily determined by the optical 
thickness, single scattering albedo and asymmetry parameter of the ice crystals. The detailed shape of the phase function is of 
second order importance. 

The optical properties of single ice particles can be calculated using GO or other methods (Yang et al., 2015). Real ice clouds 
obviously consist of numerous ice crystals that have certain distributions of sizes and shapes. Recipes for integrating optical 


properties over such distributions are given by Baum et al. (2005). 
3.2 Dependence of optical properties on crystal shape 


The optical properties of ice depend on wavelength and ice size, shape and orientation. Here, we limit the discussion to 
randomly orientated crystals. Although horizontally oriented crystals are often present in clouds, it is estimated that their 
numbers are small and their influence on the integrated cloud properties is generally minimal (Bréon and Dubrulle, 2004; 
Zhou et al., 2012). Furthermore, essentially all radiative transfer and remote sensing applications are practically limited to 
randomly oriented particles. 

The extinction efficiency mainly depends on the particle size parameter. It is commonly assumed that errors due to the GO 
approximations (1.e., Q. = 2) are sufficiently small for most applications for particles with a size parameter greater than 100 
(Baran, 2009; Bi et al., 2014). Following Bryant and Latimer (1969), here the size parameter x at wavelength \ is defined as 
an V 
A Ap’ 


where V is the volume of the particle and m,. is the real part of the refractive index of ice. This definition of size parameter 


x= (m,-1) (10) 


depends on the volume-to-area ratio of the crystals. Confusingly, various other definitions of the size parameter that depend on, 
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e.g., maximum diameters or the radius of an equivalent-area sphere, are used in the available literature (e.g., Wyser and Yang, 
1998; Zakharova and Mishchenko, 2000; Um et al., 2015). Bryant and Latimer (1969) showed that Q, is a function of y with a 
minimal dependency on shape and refractive index when the size parameter is defined as in Eq. 10. Occasionally ice extinction 
efficiencies are shown to be seemingly substantially dependent on shape, but at least part of the shape-dependency of the 
ice extinction efficiency reported in the literature is caused by using a definition of size parameter other than Eq. 10 (e.g., 
in Zakharova and Mishchenko, 2000; Yang et al., 2013; Liu et al., 2014; Um and McFarquhar, 2015). For example, extinction 
efficiencies presented as a function of maximum dimensions for different shapes may show functional differences, but the 
differences will be small if extinction efficiencies are instead presented as a function of the size parameter defined as by Eq. 10 


(cf. Wyser and Yang, 1998). 


The influence of particle shape on single scattering albedo is also limited (Wyser and Yang, 1998; Key et al., 2002; van Diedenhoven et al 


2014a) since it mainly depends on the absorption size parameter as defined by van Diedenhoven et al. (2014), i.e., 


mV 
abs =~) 11 
Xab XA, (11) 


where ™m, is the imaginary part of the refractive index. For a distribution of ice crystals, the integrated single scattering albedo 


at a given wavelength is mainly determined by the particle size distribution effective radius defined as 


_ 3 Viot 
7 4 Atot’ 


Rest (12) 


where Vio¢ and Ajo¢ are the total volume and projected area of all particles in the cloud volume (Baum et al., 2005). This domi- 
nant influence of size, rather than shape, on the single scattering albedo make measurements that depend on the single scattering 
albedo, i.e. at shortwave infrared wavelengths, very suitable in retrieving particle effective radius (Nakajima and King, 1990). 

In turn, the ice crystal scattering phase matrix is largely determined by particle shape, rather than size, especially at solar 
wavelengths where ice absorption is minimal. In fact, in the GO approximation and at non-absorbing wavelengths, the ray 
tracing part of the phase matrix is scale invariant by definition. Scattering phase matrices for crystals with size parameters 
smaller than about 100 do feature substantial size dependencies that are not accurately modeled by classical GO but these 
can be approximated by corrections to GO and/or by other techniques (Yang et al., 2015). However, most cross-sectional area 
is generally contributed by natural crystals larger than this limit when evaluated at solar wavelengths (Baran, 2009). As can 
be inferred from Eq. 8, the phase function is also affected by the single scattering albedo for non-absorbing wavelengths. 
Confusingly, substantial size-dependencies of scattering phase matrices for large particles at non-absorbing wavelengths are 
frequently presented in publications (e.g., Key et al., 2002; Baum et al., 2011, 2014). However, these variations of the phase 
matrices with size mostly stem from the assumed crystal geometries (e.g., the aspect ratios of ice crystal components) that 
depend on size. Although geometries might be expected to change with size for natural ice crystals, it is important to separate 
the dependence of optical properties on size from their dependence on crystal shape when considering the remote sensing of 
crystal shape. The strong dependence on shape of ice phase matrix elements at solar wavelengths and the minimal dependence 
on size are exploited by the methods to remotely sense information about ice shape using signatures in the scattered light, as 


discussed in section 4. 
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Figure 3. Asymmetry parameter of hexagonal prisms at a wavelength of 865 nm as a function of aspect ratio (x-axis) and 


roughness parameter (y-axis). Figure reproduced from van Diedenhoven et al. (2012a) 


Thus, out of the fundamental optical properties, only scattering phase matrices are substantially dependent on shape. For 
the purpose of remote sensing of ice crystal shapes, it is important to determine to which shape characteristics the phase 
matrix is particularly sensitive. Early applications of ray tracing calculations already determined that the phase matrix of single 
hexagonal prisms show a strong dependence on their aspect ratio (e.g, Macke et al., 1996b). For example, the dependence 
of the asymmetry parameter at a wavelength of 865 nm on the aspect ratio hexagonal prisms is shown in Fig. 3. This figure 
also depicts the dependence of the asymmetry parameter on crystal surface roughness, as discussed later. As seen in Fig. 3, 
the asymmetry parameter is lowest for crystals with an aspect ratio near unity. Also the increase of asymmetry paramneter 
with aspect ratio is rather symmetric for plates and columns, implying that asymmetry parameter is largely determined by the 
aspect ratio ~<; as defined by Eq. 2 (van Diedenhoven et al., 2016a). Asymmetry parameters increase with aspect ratio owing 
to the increase of parallel surface areas leading to greater probability of light passing through the particle with low orders of 
refraction + reflection and a minimal change of direction (Yang and Fu, 2009). For the same reason, the increase of asymmetry 
parameter with decreasing aspect ratio a <, is slightly weaker for columns than for plates owing to the larger parallel surfaces 
of plates (cf. Macke et al., 1996b; Yang and Fu, 2009). Examples of phase functions and the degree of linear polarization for 
hexagonal ice prism with varying aspect ratios and for various complex ice habits are shown in Figs. 4 and 5, respectively. 
The applied size distributions are described by van Diedenhoven et al. (2012a). From right panels in Figs. 4 and 5, one is 
tempted to conclude that ice optical properties strongly depend on the polycrystalline structure or habit. However, on closer 


inspection, the variation of phase function between different habits is mostly caused by the variation of the aspect ratio of 
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Figure 4. Phase functions of hexagonal columns with varying aspect ratios (AR, top left) and roughness parameters (lower left) 
and phase functions of complex crystals (Yang et al., 2015) with smooth (upper right) and severely roughneded (o = 0.5, lower 


right) surfaces. 


their components. For example, aggregates of columns have components with an average aspect ratio of 1.5 (Fu, 2007) and 
their phase matrix elements show resemblances with compact hexagonal columns (see Figs. 4 and 5), while bullet rosettes 
have thinner components and their phase matrices are similar to those of thin columns. Moreover, as was recognized by many 
authors (e.g., Iaquinta et al., 1995; Um and McFarquhar, 2007, 2009; Fu, 2007; Baran, 2009; van Diedenhoven et al., 2012a, 
2014a), the phase function of an ice crystal consisting of multiple, identical hexagonal components largely resembles the phase 
function of the individual single components. Comparisons of phase functions of complex crystals and their components are 
shown previously for bullet rosettes and their aggregates (Iaquinta et al., 1995; Um and McFarquhar, 2007) and aggregates of 
plates (Baran, 2009; Um and McFarquhar, 2009) and columns (Baran, 2009). For various crystal habits, Fu (2007) showed 


that the asymmetry parameters of the complex crystals are approximated to within about 0.01 by the asymmetry parameters 
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Figure 5. Same as Fig. 4 but for the degree of linear polarization. 


of their individual components. Um and McFarquhar (2007) showed that the asymmetry parameter of individual bullets at 
visible wavelengths are about 0.01 larger than rosettes with 6 arms, and about 0.02 larger than an aggregate of bullet rosettes 
consisting of 18 components. For aggregates of plates, Um and McFarquhar (2009) showed that the asymmetry parameter 
at visible wavelengths decreases with roughly 0.01 for every 5 plates attached to an aggregate, although the decrease was 
found to be non-monotonic. Furthermore, they found that differences between asymmetry parameters of aggregates of plates 
and those of the component plates increase as the distances between the centres of mass of the plates within the aggregates 
decrease. Essentially, increasing the compactness of aggregate of plates leads to stacking of the plates and finally to morphing 
into a thicker single plate and hence to a lower asymmetry parameter. To put the sensitivity of the asymmetry parameter to 
aggregation into perspective, a small change in the aspect ratio of a single plate, from 0.6 to 0.5, increases the asymmetry 
parameter at a wavelength of 550 nm by 0.01, i.e., from 0.79 to 0.80. (Here the parameterization of van Diedenhoven et al. 


(2014a) is used to calculate the asymmetry parameter.) Thus, although the larger-scale structure of such complex particles 
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influences their single-scattering properties, the influence of natural variations in aspect ratios of crystal components may be 
expected to be greater. 

As discussed in section 2.3, aside from general shape of ice crystals and the aspect ratios of their components, microscale 
structures on the crystal surfaces are of importance for optical properties. In order to compute the optical properties of a given 
ice crystal, an idealized definition of its facets is needed. However, the computational effort sharply increases as the number of 
facets required to define a particle with microscale surface roughness increases (Liu et al., 2013). Moreover, these complex or 
rough particles are determined by stochastic processes and all particles in a cloud volume essentially have unique complexity 
or roughness structures. To calculate the optical properties of such collections of particles, basically the optical properties of 
all many random members have to be determined and then appropriately averaged. Since this is impractical, stochastic ap- 
proaches to approximate the effects of particle roughening have been developed. A commonly used method is based on the 
work of Peltoniemi et al. (1989), Macke et al. (1996b) and Yang and Liou (1998) and takes crystal roughness of ice crystals 
into account in a statistical manner by perturbing the normal of the crystal facet surface from its nominal orientation by an angle 
that, for each interaction with a ray, is varied randomly with a given distribution. This was found to be an effective and accurate 
treatment of the effects of surface roughness on ice crystal scattering properties (Liu et al., 2013). Practically, codes that use 
this approach (Macke et al., 1996b; Yang and Liou, 1998; Shcherbakov et al., 2006) vary in the assumed distribution of facet 
tilt angles (Neshyba et al., 2013; Geogdzhayev and van Diedenhoven, 2016). For instance, the code by Macke et al. (1996b) 
varies the tilt angles randomly with uniform distribution between 0° and 6 - 90°, where 6 is referred to as the roughness pa- 
rameter. Yang and Liou (1998) and Shcherbakov et al. (2006) vary tilt angles according to Gaussian and Weibull distributions, 
respectively, where the mean tilt angle is determined by roughness parameter 0. However, Geogdzhayev and van Diedenhoven 
(2016) and Neshyba et al. (2013) showed that roughness parameters with these different definitions, for a given value, lead to 
rather similar phase functions and asymmetry parameters (i.e. 6 + 7). However, these roughness parameterizations currently 
need to be considered as merely effective methods to simulate the effects of roughness structures on the scattering matrices, 
since they are difficult to relate to physical characteristics of the structures observed on real ice crystals, although attempts 
have been made (Neshyba et al., 2013). Furthermore, Liu et al. (2014) have shown that the simulated effects of microscale 
roughness features and larger scale geometric irregularities of ice crystals are similar. Indeed, irregular crystals such as Koch 
fractals (Macke et al., 1996b) and Voronoi particles (Letu et al., 2016) have phase matrices similar to rough hexagonal prisms. 
Increasing the number of impurities within ice crystals also has a similar effect (Macke et al., 1996a; Hong and Minnis, 2015; 
Panetta et al., 2016). Thus, the roughness parameterization can be interpreted as representing the randomization of reflection 
and refraction on and in the particles caused by any of such deviations from pristine, solid, smooth crystals. Since the roughness 
parameter reflects more than just roughness structures, some authors prefer the term "distortion parameter" instead. Using the 
definition of roughness of Macke et al. (1996b), Fig. 3 shows the dependence of the asymmetry parameter on roughness and 
aspect ratio, while the bottom left panel of Figs. 4 and 5 show the variation of phase functions and degree of linear polariza- 
tion, respectively, on roughness for columns. Also, the bottom right panels of Figs. 4 and 5 show the phase functions and the 
degree of linear polarization, respectively, of severely roughened complex particles (o = 0.5). Angular features in the phase 


functions and degree of linear polarization functions are smoothed out as roughness increases. Most notable, the halos are 
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diminished as roughness increases and asymmetry parameters systematically decrease with roughness (cf. Yang et al., 2008a; 
van Diedenhoven, 2014; van Diedenhoven et al., 2014a; Geogdzhayev and van Diedenhoven, 2016). 

Another property of ice crystals that can have substantial impact of scattering properties is the presence of hollow crystal 
endings, or cavities (Schmitt et al., 2007; Yang et al., 2008b; Smith et al., 2015). Generally, asymmetry parameters increase as 
cavity depth increases, albeit this increase becomes smaller for crystals with aspect ratios increasingly deviating from unity 
(Schmitt et al., 2007). As can be seen in Fig. 4 (right panels), the hollow structures increase the scattering phase function values 
at scattering angles smaller than 20 degrees, but decrease values at side- and backscattering angles of 60 to 180 degrees (cf. 
Schmitt et al., 2007; Yang et al., 2008b). One way to interpret the observed influence of hollowness on columns is that, since 
the walls of the hollow parts of the crystals are thin, these essentially resemble columns with high aspect ratios, leading to an 
increased asymmetry parameter (van Diedenhoven et al., 2012a). This interpretation is also supported by the resemblance of 
phase functions and degree of linear polarization of hollow particles to those of columns with large aspect ratios seen in Figs. 
4 and 5. 

In summary, the most important features of ice particle shape influencing their optical properties are the aspect ratios of 
the crystals or the components of complex crystals as well as the microscale surface roughness, distortion or hollowness. The 
larger scale polycrystaline structure of the ice crystal is of lesser importance. Thus, the different characteristics of ice particle 


shape can be grouped in order of importance for the scattering properties as follows: 
1. Aspect ratios of simple crystals or components of complex crystals; 
2. Microscale structure, surface roughness, distortion, impurities or cavities; 
3. Polycrystalline structure or habit. 


These crystal shape characteristics are relevant targets for the remote sensing techniques described in the next section. 

As discussed in section 2, myriads of shapes are found in real cloud volumes. Moreover, the aspect ratios and roughness levels 
may be expected to vary per particle (Um et al., 2015; Fridlind et al., 2016; Schnaiter et al., 2016). However, most publications 
about optical properties of ice crystals and remote sensing of ice crystal shapes consider calculations for individual crystals 
with a particular shape rather than mixtures of crystals with large variations in shapes, aspect ratios and roughness levels. 
Various habit mixtures have been presented (e.g., Baum et al., 2005, 2011; Baran and C.-Labonnote, 2007) but these usually 
contain just a handful of ice shapes. Furthermore, remote sensing applications generally aim to find a single shape or simple 
mixture that fits the measurements best. The question arises how to interpret remote sensing results that conclude one crystal 
shape or a limited mixture to be consistent with the measurements. In other words, which part of the ensemble is represented by 
the single inferred crystal shape? In order to estimate the asymmetry parameter of an ensemble of hexagonal ice crystals with 
a distribution of aspect ratios, Fu (2007) defined an ensemble-average aspect ratio as the average aspect ratio weighted by the 
orientation-averaged projected area of the crystals in the distribution. van Diedenhoven et al. (2016a) showed that the definition 
of aspect ratio limiting the values below unity for both plates and columns (i.e., Eq. 2) needs to be used in order for an ensemble- 


average aspect ratio to adequately represent the optical properties of the ensemble. Furthermore, they showed that even for 


17 


10 


15 


20 


25 


30 


ensembles containing both columns and plates, the ensemble asymmetry parameters are generally represented consistently 
by a single crystal with an aspect ratio equal to the ensemble average, especially if geometrical averaging is used. Tests 
on mixtures of plates and columns yielded root-mean-squared differences between ensemble-average asymmetry parameters 
and those calculated from ensemble-average aspect ratios were 0.006. Furthermore, van Diedenhoven et al. (2016a) showed 
that effective asymmetry parameters based on arithmetic averages of roughness parameters are also largely consistent with 
ensemble-average asymmetry parameters. As discussed by Fu (2007) and van Diedenhoven et al. (2016a), these conclusions 
are likely qualitatively applicable to ensembles of more complex structures such as aggregates of columns, aggregates of plates 
and bullet rosettes with a range of aspect ratios of their components, as well as to internal mixtures of plate-like and column-like 
components, such as asymmetric bullet rosettes, aggregates of varying plates and/or columns or plate-capped columns. Thus, 
the conclusions by Fu (2007) and van Diedenhoven et al. (2016a) suggest that shapes inferred from remote sensing applications 
represent an area-weighted ensemble average of the ice crystals in the observed cloud volume. Further discussion on this is 
included in section 4.2.4. 

In summary, it is practically impossible to consider all geometries of ice crystals formed in natural ice clouds for remote 
sensing applications. Any ice model used for remote sensing purposes will be a highly idealized simplification of natural ice 
crystals. Often, crystal shapes that are found to match observations are presented merely as "radiative equivalent effective 
shapes" (e.g., McFarlane and Marchand, 2008; Cole et al., 2013), and no direct relation between the inferred shape and the ice 
crystal shapes in the observed clouds is claimed. However, using idealized models, it is still possible to infer useful information 
about shape characteristics of natural ice crystals in clouds. Here, a famous quote from the late George Box is fitting, namely 
that "essentially, all models are wrong, but some are useful" (Box and Draper, 1987). In order to infer information about 
the variation of ice crystal shape in ice clouds, a set of ice models is needed that systematically varies the physical and 
scattering properties and is able to fit the full range of available remote sensing observations. George Box provides further 
inspiration by stating that "since all models are wrong, the scientist cannot obtain a ‘correct’ one by excessive elaboration. On 
the contrary, following William of Occam [she/]he should seek an economical description of natural phenomena" (Box, 1976). 
Following such wisdom, one can conclude that for remote sensing it may be advised to focus on the variation of the ice shape 
characteristics that are mostly determining the radiative properties as listed above, i.e., component aspect ratios and microscale 
structure, and less on crystal habits. Note that shape characteristics such as aspect ratios, roughness levels and cavity depths 
are quantitative parameters that suit systematical variation in contrast to *habit’. Focusing on such shape parameters provides 
a "economical description" and reduces the remote sensing problem substantially to quantifying these parameters without 
considering the virtually infinite number of possible large-scale shapes. Based on the equivalence of scattering properties 
of complex crystals and their components, as discussed above, such inferred parameters, e.g., aspect ratios and roughness 
parameters, represent the averaged properties of the components of the crystals in the observed ice clouds. Given the variation 
of ice crystals observed and the myriad of conditions in which ice crystals form and evolve, we can assume that values of aspect 
ratios and roughness levels occur continuously over rather large ranges, which need to be spanned by the models considered 


for remote sensing applications to infer ice crystal shapes. 
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4 Remote sensing of ice crystal shapes 


Essentially two different observational approaches have particular potential for inferring information about particle shape, 
namely active lidar observations and passive multi-directional measurements of total and/or polarized reflectances. Here, these 


approaches are discussed separately in sections 4.1 and 4.2. 
4.1 Lidar measurements 


Lidars (i.e., laser radars) probe clouds by emitting laser beam pulses and measuring the power of the returned signals (Weitkamp, 
2005). By measuring the delay between the emitted and detected beams, the distance between the lidar and the cloud volume 
on which the beam was scattered back can be determined, resulting in height-resolved information. When the emitted beam is 
polarized and the polarization state of the detected signal is measured, the depolarization of the beam caused by backscattering 
on particles can be measured. The altitude-dependent linear depolarization ratio 6; of the returned signal is usually defined as 
(Schotland et al., 1971) 


_ B1(z) 
By (z)’ 


where ( is the backscattering cross sections in the planes of polarization perpendicular (L) and parallel (||) to the laser’s 


61(2) (13) 


reference plane. In Eq. 13 it is assumed that atmospheric extinction is independent of the polarization state of the propagating 
radiation. Other definitions for lidar depolarization are discussed by Gimmestad (2008). 

Since lidar signals are mainly determined by singly scattering light, the depolarization of a lidar signal returned from a cloud 
layer can be simulated in a straight forward manner from the assumed scattering phase matrix, namely by 


_ P11(180°) — P22(180°) 
~ Py (180°) + Po2(180°) ° 


61(2) (14) 


Lorenz-Mie theory applied to spherical particles shows that no depolarization is caused by single back-scattering on purely 
spherical cloud droplets, although multiple scattering on liquid cloud drops can lead to some depolarization of the lidar signal 
(Sassen and Petrilla, 1986; Hu et al., 2001). However, single back-scattering of the laser light on non-spherical particles such 
as spheroids or hexagonal prisms leads to substantial depolarization of the signal that is largely dependent on particle shape. 
The first ground-based lidar depolarization measurements of clouds were reported by Schotland et al. (1971) and many po- 
larized ground-based lidars are deployed world-wide today. The satellite-based Cloud-Aerosol Lidar with Orthogonal Polariza- 
tion (CALIOP) on the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) platform (Winker et al., 
2007), launched in 2006, is particularly relevant for ice cloud studies as it is providing global statistics and observes cirrus 
otherwise obscured by lower lying liquid clouds when viewed from the surface. Similar statistics are expected from the Atmo- 
spheric Lidar (ATLID) on the Earthcare platform to be launched in 2018 (Lefebvre et al., 2016). Furthermore, the cloud physics 
lidar (CPL) has been providing a wealth of data since the year 2000 while mounted on the high altitude NASA ER-2 aircraft 
during many field campaigns targeting clouds (Yorks et al., 2011). General conclusions about the depolarization measured by 


such lidars is that the depolarization ratio in ice clouds typically range from 0.2 to 0.6, and generally increases with increasing 
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cloud height or decreasing temperature (Sassen and Benson, 2001; Reichardt et al., 2002; Noel et al., 2004; Sassen and Zhu, 
2009; Yorks et al., 2011; Baum et al., 2011; Sassen et al., 2012). Furthermore, CALIOP statistics presented by Sassen et al. 
(2012) show a globally, vertically and yearly averaged ice cloud depolarization ratio of 0.37 and generally lower values at 
higher latitudes. In addition, Martins et al. (2011) found no significant variations in depolarization ratios for ice clouds in rela- 
tion to updraft strength or horizontal windspeed. Note that the depolarization ratios quoted here reflect the values measured at 
an off-nadir angle, since horizontally orientated ice crystals strongly decrease depolarization measured at direct nadir owing to 
specular reflection on the crystal surfaces (Del Guasta et al., 2006; Zhou et al., 2012; Sassen et al., 2012). 

Although the variation of depolarization ratios in ice clouds at different altitudes, temperatures and latitudes are fairly well 
documented, a definite quantitative interpretation of the values in terms of ice crystals shape variations is lacking. One in- 
terpretation is offered by Noel et al. (2002) who showed that the simulated depolarization ratio of single hexagonal prisms 
depends on their aspect ratio. As also shown by Yang and Fu (2009), the simulated depolarization ratio of hexagonal crystals 
with smooth surfaces is about 0.2 for thin plates, oscillating between 0.3 and 0.4 for thicker plates and about 0.4—0.6 for com- 
pact particles and columns. Based on this relation between depolarization ratio and aspect ratio, Noel et al. (2004) concluded 
that the depolarization ratios observed by the CPL during the CRYSTAL-FACE campaign conducted near Florida in 2002 
are consistent with a dominance of compact or ’irregular’ crystals and an increasing contribution by plate-like crystals with 
increasing cloud top temperature. Measurements on ice crystals grown in the laboratory presented by Amsler et al. (2009) and 
Abdelmonem et al. (2011) seem to support this quantitative interpretation of depolarization ratios, although modeled depo- 
larization ratios were found to be generally larger than the measured ones. However, as discussed in section 2, most crystals 
observed in situ have more complex, aggregated shapes than simple hexagonal prisms. Although large scale complexity gener- 
ally has a weaker effect on the scattering phase matrix than aspect ratio, as discussed in section 3.2, it is unclear whether or how 
this complexity affects the depolarization ratios specifically (Reichardt et al., 2008). A fundamental problem for simulating li- 
dar backscattering and depolarization on ice crystals is that the phase matrix computed based on geometric optics exhibits a 
singularity at the exact backscattering direction for randomly oriented crystals (Borovoi et al., 2005, 2014). Usually values for 
angles close to 180 degrees are used instead or they are extrapolated towards 180 degrees (e.g., Smith et al., 2016). Further- 
more, the dependency of depolarization on aspect ratio presented previously by Noel et al. (2004) and others are based on solid 
crystals with smooth surfaces. Although calculations presented by Smith et al. (2016) suggest that the effects of crystal hol- 
lowness on lidar depolarization is generally small, reported estimates of the effects of surface roughness on lidar depolarization 
are inconclusive. Mostly, roughness or distortion of ice crystals is shown to lead to a decrease of depolarization (Del Guasta, 
2001; Baum et al., 2010, 2011), although the opposite trend has been shown as well (Smith et al., 2016; Konoshonkin et al., 
2016). Furthermore, rounded crystal shapes caused by crystal sublimation are sometimes associated with a decrease of de- 
polarization (e.g., Martins et al., 2011), although increasing depolarization from sublimation is suggested by Schnaiter et al. 
(2012). Also, low depolarization in ice clouds are sometimes interpreted as caused by compact or ’quasi-spherical’ particles 
(e.g., Choi et al., 2010). However, these associations may not be supported by observations and calculations of optical prop- 
erties. For example, Mishchenko and Sassen (1998) showed that depolarization ratios of spheroids and deformed spheroids 


with aspect ratios close to unity are comparable with, or larger than, the values of hexagonal plates, even if such shapes can 
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be considered quasi-spherical. Sassen (1977) reported depolarization ratios above 0.5 for frozen rainwater drops with a reg- 
ular spherical or spheroidal appearance. In contrast, recent work presented by Jarvinen et al. (2016) and Schmitt et al. (2016) 
shows that sublimation of ice crystals can lead to smooth frozen droplets that have optical properties similar to spheres includ- 
ing near-zero depolarization ratios. Other factors complicating the interpretation of lidar depolarization values are the presence 
of oriented ice crystals (Del Guasta et al., 2006; Zhou et al., 2012) and the co-existence of ice crystals and liquid drops in 
mixed-phase clouds (Bourdages et al., 2009; van Diedenhoven et al., 2011) 

In conclusion, systematic variations of lidar depolarization in clouds are observed, but their interpretation is currently incon- 
clusive. As put by Sassen et al. (2012), "these [variations in depolarization ratios] must reflect the different ice particle shapes 
that depend on the cirrus cloud formation mechanism (e.g., convective-anvil, orographic, synoptic, etc.), or more specifically on 
the basic cloud-particle forming aerosol available for crystal formation and the microphysical effects of typical cloud updraft 
velocities." However, currently, interpretations of lidar depolarizations in ice clouds are more qualitative than quantitative. To 
quantitatively relate the measured depolarization ratios to microphysical properties of the ice crystals, more research is needed, 
including 1) further development of methods for calculating exact backscattering properties of complex ice crystals and crys- 
tals with randomly roughened surfaces; 2) investigations on the effects of crystal complexity on lidar depolarization; 3) further 


statistical studies on measured lidar depolarization ratios with co-located in situ ice crystal shape observations. 
4.2 Multi-angular measurements 


As shown in Figs. 4 and 5, the shapes of the scattering phase function and degree of polarization of ice crystals for visible 
wavelengths depend on the crystal shape. These angular variations in scattering phase function and degree of polarization lead 
to angular features in the total and polarized reflectances of ice clouds that can be observed by remote sensing instruments that 
make observations at multiple viewing angles per instrument footprint. Here, we first sumarize the techniques used to infer 
ice crystal shape from multi-angular total reflectances (section 4.2.1) and multi-angular polarized reflectances (section 4.2.2). 
Then some discussion about data selection and availability is offered in section 4.2.3 before summarizing results in section 


4.2.4. 
4.2.1 Multi-angular total reflectances 


The measured reflectance of a cloud layer can be defined as 


ve A 
R10, Ad) = “ite (15) 


where /49 and yu are the cosine of the solar zenith angle and viewing angle, respectively, Fi, is the solar irradiance, and I is the 
radiance measured by an instrument. Since this definition refers to the reflectance of light irregardless the polarization state, it 
is commonly referred to as "total reflectance". 

Ice crystal shape can be inferred from total reflectances because angular features in the single scattering phase function of 
ice crystals, as seen in Fig. 4 for example, are preserved in the directional measurements. The directional reflectance of singly 


scattered light is determined by the phase function (van de Hulst, 1957). Because of the large extinction in cloud tops, however, 
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singly scattered light emerging from the cloud is only a small fraction of the total reflection, adding a weak signal of direc- 
tionality to the measured total reflectances. As discussed in section 3, the magnitude of the total reflectance of an ice cloud is 
mostly determined by its optical thickness and the asymmetry parameter of the scattering phase function (Coakley and Chylek, 
1975). One might expect that multiply scattered light emerging from a cloud is fully isotropic. However, as explained by 
Zhang et al. (2009), at any scattering event, the highly peaked phase functions of ice crystals lead to a large portion of the rays 
scattered forward without a substantial change in direction. Thus, a likely path for light rays is one where one or more forward 
scattering events are followed by a single side-ward scattering event, which is then subsequently followed by one or more 
forward scattering events again before light emerges from the cloud and is measured by an instrument. Although such paths 
describe multiply scattered light, the angular features in the scattering matrix of ice crystals are also preserved in the directional 
reflectance measurements for these paths, similarly as for the singly scattered light. This effect enhances the signal to noise 
ratio and the information content of multi-directional total reflectance measurements for the use of retrieving ice crystals shape. 

For a plane-parallel homogeneous cloud layer, its visible reflectance at a particular viewing geometry is fully determined 
by the scattering phase function and the optical thickness. Thus, when multi-directional measurements are available for a 
cloud layer, an optical thickness can be retrieved from each angular measurement assuming a particular phase function. If the 
assumed phase function is correct then the optical thickness retrieved at each angle would be exactly the same, i.e., the correct 
value. However, an incorrect phase function assumption leads to an angular dependence of the retrieved optical thickness. 
This angular variation of retrieved optical thickness, or spherical albedos derived from them (e.g., Doutriaux-Boucher et al., 
2000), is the basis of methods to retrieve ice shapes from multi-angular total reflectances (e.g., Doutriaux-Boucher et al., 
2000; McFarlane and Marchand, 2008). Generally, the angular variation of retrieved optical thicknesses or spherical albedos 
is determined for a selection of assumed crystal shapes and crystals leading to the smallest angular variations are presented as 
the most consistent with real ice crystals. 

Examples of instruments that make measurements of total reflectances at multiple viewing angles per footprint are NASA’s 
Multi-angle Imaging SpectroRadiometer (MISR, Diner et al., 2002) and ESA’s Along-Track Scanning Radiometer 2 (ATSR-2) 
and the Advanced ATSR (AATSR) (Sayer et al., 2011). The ATSR-2 and AATSR instruments have nadir pointing cameras as 
well as one forward pointing camera observing each footprint. MISR has 9 cameras at nadir, forward and aft pointing angles. 
Also multi-angle polarimeters as discussed in section 4.2.2 measure total reflectances in addition to the polarized reflectances. 
Overlapping footprints of geostationary satellite imagers have also been used to create dual-view observations (Chepfer et al., 
2002). Furthermore, global statistics of single-view imager data can be used to infer the angular variation of cloud reflectances 


(Zhang et al., 2009; Wang et al., 2014). 
4.2.2 Multi-angular polarized reflectances 


Similar to the definition of total reflectances (Eq. 15), polarized reflectances R,, are generally defined as 


VF 
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The dependency of Q and U on po, ps and A¢ is omitted in Eq. 16 for clarity. Note that both Q and U are signed according 
to the orientation of the measured electromagnetic wave relative to a reference plane. This information about polarization 
direction is lost by the definition of Eq. 16, but can be analyzed separately (e.g., Sun et al., 2015) or a sign can be added to the 
polarized reflectances accordingly, as described by C.-Labonnote et al. (2001). Normalizations different from the one used in 
Eq. 16 are sometimes applied too (e.g., C.-Labonnote et al., 2001). 

Ice crystals shape characteristics can be inferred from the polarized reflectances by matching the measured angular variation 
of polarized reflectances with values simulated using a radiative transfer model that includes multiple scattering and polar- 
ization. Since multiple scattering depolarizes light, polarized reflectances are dominated by singly scattered light and have 
pronounced angular dependencies that are determined by the polarization properties of the ice crystals, as shown in Fig. 5. Fur- 
thermore, for cloud optical thickness values above about 2—5 (depending on the asymmetry parameter of the phase functions), 
polarized reflectances do not depend on cloud optical thickness anymore (van Diedenhoven et al., 2012a), which is convenient 
when retrieving information about ice crystals shape in various thick cloud types. However, measuring the polarization state of 
reflected light is challenging, as per footprint and per viewing angle it requires a minimal set of three simultaneous observations 
at three well-determined polarization angles (Tyo et al., 2006). 

Relevant multi-angular polarized reflectances from space were measured by the POLarization and Directionality of the 
Earth’s Reflectances (POLDER) instruments (Deschamps et al., 1994; Fougnie et al., 2007). The POLDER imager acquired 
overlapping images at different orbital locations that yield observations for each given footprint at a maximum of 16 viewing 
angles (Fougnie et al., 2007). For the POLDER wavelength bands at 490, 670, and 865 nm, three acquisitions are performed 
through a polarizer oriented at -60°, 0°, and 60° relative to a given reference from which the Stokes parameters [I,@Q,U] are 
retrieved. In addition to an airborne version of the POLDER instrument (Chepfer et al., 1998), three satellite-based versions 
were deployed. The first two (POLDER-1 and POLDER-2) were mounted on the Advanced Earth Observing Satellite (ADEOS- 
1 and ADEOS-2, respectively), but both only lasted for about 6 months owing to platform failures. The third, mounted on the 
CNES/Myriade Polarization & Anisotropy of Reflectance for Atmospheric Sciences coupled with Observations from a Lidar 
(PARASOL) microsatellite, was highly successful and lasted for nearly 9 years, from 2004 to 2013. POLDER/PARASOL 
was part of NASA’s afternoon train (A-train) constellation until 2010, allowing its data to be combined with that of Cloudsat, 
CALIPSO, MODIS and other A-train instruments. No multi-directional polarimeter other than POLDER has flown in space 
to date. Unfortunately, the launch of the Glory satellite carrying the Aerosol Polarimetry Sensor (APS, Mishchenko et al., 
2007) failed in March 2011. Several airborne multi-angle polarimeters are regularly deployed during field campaigns, often on 
high altitude aircraft capable of flying sufficiently above ice clouds. The Research Scanning Polarimeter (RSP, Cairns et al., 
2003) is an airborne version of the APS and provides simultaneous measurements of total and polarized reflectances in 9 
spectral bands from the visible to the shortwave infrared. The RSP scans along track, providing measurements of each pixel 
at 152 different viewing angles at 0.8 degrees intervals. Other deployed airborne polarimeters include the Airborne Multi- 
angle Spectro Polarimetric Imager (AirMSPI, Diner et al., 2013) and Airborne Spectropolarimeter for Planetary EXploration 
(SPEX, Rietjens et al., 2015). 


23 


20 


25 


30 


4.2.3 Data selection 


For the remote sensing of ice crystal shape, employing relevant data selection criteria is crucial. Obviously, one impor- 
tant selection criteria is that ice clouds are selected without interference of liquid clouds. For studies on POLDER data, 
this is generally efficiently achieved by testing for rainbow features in polarized reflectances. However, as noted by, e.g, 
C.-Labonnote et al. (2001) and Cole et al. (2014) this generally filters out clouds with low optical thicknesses (below optical 
thickness 1-5). This inadvertent selection of thick clouds is often not noted by studies using POLDER data, although evi- 
dently statistics of POLDER polarized reflectances generally show saturated reflectances, which only occur at about optical 
depths of about 4—5 (Chepfer et al., 2001; van Diedenhoven et al., 2013). Since cirrus typically has an optical thickness below 4 
(Rossow and Schiffer, 1999), this selection has the important implication that that most cirrus are not included in the studies on 
ice crystals shapes using POLDER data. Alternatively, collocated MODIS thermodynamical phase determination can be used 
to select ice-topped clouds. However, this cloud phase retrieval is shown to be less reliable for mixed-phase and multi-layered 
clouds (e.g., Riedi et al., 2010). Interference of liquid drops to the multi-angular total or polarized reflectances is expected to 
lead to spurious angular dependence of the measurements and thus to biases the retrievals of ice shape towards more pristine 
particles. 

Since the simulations of ice optical properties and radiative transfer calculations generally assume random orientation of 
the ice crystals, Chepfer et al. (2001) and Sun et al. (2006) concluded that the POLDER data needs to be screened for the 
presence of oriented ice crystals. Specular reflection on oriented crystals is highly polarizing and leads to sharp angular peaks 
in the polarized reflectances at the specular reflection angles (Chepfer et al., 1998; Bréon and Dubrulle, 2004). This peak can 
be used to screen the data for particle orientation, but it requires measurements at and around the scattering angle associated 
with specular reflection for a given solar angle. As noted by Sun et al. (2006), this requirement excludes most data, especially 
for instruments with relatively low angular resolution such as POLDER. For example, Sun et al. (2006) noted that the stringent 
data selection criteria used in their study on POLDER data lead to only 0.37% of the 5-month POLDER dataset passing these 
criteria. However, analyzes of global POLDER data suggests that that the relative number of oriented crystals is generally low 
and that the typical effective contribution to particle area from oriented plates in clouds is smaller than 1% (Bréon and Dubrulle, 
2004). Thus, the total and polarized reflectances from ice clouds outside of the specular reflection geometry are dominated by 
randomly oriented crystals, which implies that the screening for random orientation is not needed. It may be needed to remove 
data affected by sunglint on ocean surfaces (van Diedenhoven et al., 2013), although this will only be an issue for optically thin 
clouds that are generally already excluded from POLDER measurements. 

As the crystal shape information is inferred from angular variations in total and polarized reflectances, the information con- 
tent depends on the scattering angles that are sampled by the instrument. In turn, the sampled scattering angles are dependent 
on the solar zenith angle and thus on latitude when analyzing data from polar orbiting satellites. It is important to note that 
for the statistics on polarized reflectances derived from global POLDER measurements, data from higher latitudes contribute 
more towards the small scattering angles, while tropical clouds contribute more to the larger scattering angles (Buriez et al., 


2001). This also implies that any latitudinal variations in retrieved ice model or roughness level may stem from the latitudinal 
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variation of available scattering angles and the sensitivity of the retrieval approaches to scattering angle range. Results from 
van Diedenhoven et al. (2012a, 2013) show that an angular range including samples between scattering angles of at least 120° 
and 150° is needed for ice crystal shape retrievals. However, the retrieval performance is shown to be rather insensitive to num- 
ber of angles sampled within that range or to any reasonable random error or bias on the measurements (van Diedenhoven et al., 


2012a). 
4.2.4 Overview of results 


An overview of the most notable studies on inferring ice crystal shape along with the targeted locations and the data used is 
given in Table 1. Some studies use total reflectances R, some use polarized reflectances R, and some use both. Most studies 
analyze data averaged over time and space, while some infer crystal shapes from data at the instrument’s pixel resolution. The 
table also lists a description of the dominating ice shape found by the particular studies. From this overview it is clear that 
most data are fitted best by crystals with roughness, impurities or other distortions. Some studies, especially early ones, did not 


include particle roughness or other distortions and their results should be treated with caution. 
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Reference(s) 


Region or target 


Data used 


Table 1: Overview of studies on remote sensing of ice crystal shapes 


Dominating shapes 


Notes 


Baran et al. (1998, 1999) 


Tropical anvil, cirrus 


ATSR-2 dual view 


Koch fractals 


No roughness or 


inclusions 


Chepfer et al. (1998) 


EUCREX’94 


paign, mid-lat. cirrus 


cam- 


Airborne POLDER 


Rand R, 


Pristine thin plates, 


a= 0.05-0.1 


No roughness or 


inclusions 


Doutriaux-Boucher et al. 


(2000); POLDER-1 R and | Columns with inclu- 
Global ocean & land 
C.-Labonnote et al. Rp sions, a=5 
(2000, 2001) 
Koch fractals and 
Pixel level | columns at low lati- | No roughness or 
Chepfer et al. (2001) Global ocean & land 
POLDER-1! R, tudes; Plates at high | inclusions 


latitudes 


Chepfer et al. (2002) 


Continental USA 


GOES West and East 


dual view 


Compacts, columns, 


and bullet rosettes 


No roughness or 


inclusions 


Knap et al. (2005) 


Global ocean 


POLDER-2 
Ry 


R and 


Columns with inclu- 
sions or roughness, a 


=2.5 


Knap et al. (2005) 


Tropical anvil out- 


ATSR-2 dual view 


Moderately rough- 


flow ened columns 
Pixel level | Plates and hollow | Few rough parti- 
Sun et al. (2006) Global ocean 
POLDER-1 R,, columns cles considered 
Baran and C.-Labonnote POLDER-2 R and | Rough crystals, 6 = 
Global ocean 
(2006) Rp 0.4 
McFarlane and Marchand | Southern Great | Pixel level | Rough aggregate and 
(2008) Plains MISR+MODIS bullet rosettes 
Rough plates: a=0.7 | Only 2 plate as- 
van Diedenhoven et al. POLDER- 
TWP-ICE campaign at T<-38°C; a=0.15 | pect ratios con- 
(2012b) PARASOL R, 


at T>-38°C 


sidered 
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Table 1: (continued) 


Reference(s) Region or target Data used Dominating shapes | Notes 
van Diedenhoven et al. CRYSTAL-FACE Rough compact and 
Pixel level RSP R, 
(2013) campaign Florida plate-like crystals 
POLDER- 
Cole et al. (2013) Global ocean PARASOL R | Rough habit mixture 
and R, 
Properties var 
Pixel level : 
van Diedenhoven et al. Rough plates with | with height 
TWP-ICE campaign | POLDER- 
(2014b) 6=0.4—0.7 and convective 
PARASOL R, 
strength 
Pixel level 
Rough aggregate of | Roughness varies 
Cole et al. (2014) Global ocean POLDER- 
columns with latitude 
PARASOL R, 


Wang et al. (2014) 


Global ocean & land 


MODIS over opti- 


cally thin cirrus. 


Mixture of rough and 


smooth crystals 


Differences _be- 
tween ocean and 


land 


POLDER- 
Baum et al. (2014) Global ocean Rough habit mixture 
PARASOL R, 
Data interpreted 
Pixel level 
Rough crystal mix- | as roughness 
Baran et al. (2015) Off coast of Scotland | POLDER- 
tures varying with 
PARASOL R 
humidity 
POLDER- 
Letu et al. (2016) Global ocean "Voronoi" habit 
PARASOL R 
Unphysically 
Pixel level | Rough aggregate of 
oe _ | large roughness 
Hioki et al. (2016) Global ocean POLDER- columns (only habit 
parameters for 
PARASOL R, considered) 


74% of data 
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Many of the studies listed in Table | aim to find ice crystal models that best represent globally averaged data. Some stud- 
ies compare total reflectances with modeled reflectances in order to identify crystal shapes that lead to most realistic angular 
features on average (e.g., Doutriaux-Boucher et al., 2000). Studies on global polarized reflectances usually present observa- 
tion density plots of the POLDER measurements as a function of scattering angles to which simulated polarized reflectances 
are compared (e.g., Knap et al., 2005; Baran and C.-Labonnote, 2006; Cole et al., 2013). These density plots represent the 
statistics of polarized reflectances observed at any specific scattering angle. The range of measured polarized reflectances can 
be assumed to stem from the natural variation in ice crystals shapes. It is generally concluded that the ice shapes that lead 
to simulated polarized reflectances that are closest to most observations at all of the scattering angles is a good model to 
represent natural ice clouds. However, it is important to note that almost all models included in such studies generally lead 
to simulated polarized reflectances that fall within the observation envelopes at most scattering angles. This implies that all 
such models could represent a subset of the measurements. For example, Cole et al. (2013) evaluate simulated polarized re- 
flectances for habit mixtures with smooth or rough surfaces, as well several habits with severely rough surfaces (droxtal, solid 
3D bullet rosette, hollow 3D bullet rosette, hollow column, solid column, plate, compact aggregate of columns, small spatial 
aggregate of plates, and large spatial aggregate of plates) and all simulations fall within the observation envelopes. The results 
of Baran and C.-Labonnote (2006) show that polarized reflectances simulated assuming smooth bullet rosettes and smooth 
chain-like aggregates of columns, as well as rough and distorted crystals, fall within the global statistics of measured polar- 
ized reflectances. In addition, the range of observed angular variations of total reflectances is also large in comparison to the 
variation in modeled angular features for different ice models (Doutriaux-Boucher et al., 2000; Baran and C.-Labonnote, 2006; 
Baran, 2009; Cole et al., 2013), again suggesting that many considered models may fit well to a subset of the measurements. 
Many such studies aim to select and test optical models for global retrievals of ice cloud optical thickness and effective radius 
and selecting a model that fit best to most of the globally averaged data is very much justified in this case. However, as pointed 
out by McFarlane and Marchand (2008) "using a featureless phase function will likely result in the correct scattering properties 
on average, however individual cases may have large errors." To avoid such biases, simultaneous retrievals of ice crystal shape 
and size can be employed on a pixel level (McFarlane and Marchand, 2008; van Diedenhoven et al., 2014b). For the purpose 
of collecting information about how ice crystal shape varies with, e.g., temperature, cloud type and atmospheric state, it is 
important to study the variation in ice shape and crystal roughness leading to the ranges of observed angular total and polarized 
relfectances. 

Studies that focus on the variation of ice crystals shapes are presented by, e.g, Chepfer et al. (2001), Sun et al. (2006), 
Cole et al. (2014) and Baran et al. (2015). Analyzing global POLDER-1 data, Chepfer et al. (2001) concluded that "poly- 
crystals" [i.e., Koch fractals (Macke et al., 1996b)] and hexagonal columns seem to dominant at low latitudes, whereas the 
hexagonal plates seems to occur more frequently at high latitudes. However, other than the Koch fractal, no other distorted or 
roughened crystals are included in that study. More recently, Cole et al. (2014) included 9 different shapes (Yang et al., 2015) 
and a mixture (Baum et al., 2011) with a large range of simulated roughness levels in their global study. They found that the 
aggregate of columns dominates at all latitudes and plates are the least representative of the POLDER-PARASOL polarized 


reflectances globally. Interestingly, particles with smooth surfaces (no or low roughness levels) were found to be more prevalent 
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at high latitudes, while severely rough crystals (o=0.5) were observed most frequently in the Tropics. This suggests that, in 
general, ice phase functions and polarization functions have some more angular features in the high latitudes than at low lati- 
tudes, which may be consistent with the findings of Chepfer et al. (2001) that smooth hexagonal plates fit the POLDER-1 data 
more frequently at high latitudes than elsewhere. Few roughness parameter values o > 0.5 were found by Cole et al. (2014). 
Somewhat in contrast with these results are the results of Hioki et al. (2016). Applying an algorithm based on empirical or- 
thogonal function analysis of the modeled and measured POLDER-PARASOL polarized reflectances and assuming aggregates 
of columns with varying roughness, Hioki et al. (2016) found roughness values varying with latitude, but they obtained o > 1 
for most data, which is unrealistically large. They conclude that the unexpected results indicate that the roughness retrieval is 
sensitive to an assumed particle shape, although the same aggregates of columns model with o < 0.5 is found to be matching 
the data well by Cole et al. (2014). The results by Sun et al. (2006), analyzing global POLDER-1! data, also appear to be in 
conflict with the studies by Chepfer et al. (2001) and Cole et al. (2014) and arguably every other study listed in Table 1, since 
rough particles were found to fit virtually none of the measurements and smooth plates and hollow columns were inferred from 
the data most frequently. Smooth plates and hollow columns lead to strong angular variation of polarized reflectances that are 
not often seen in the POLDER data and it may hence be surprising that Sun et al. (2006) found these habits to fit most of the 
included data. The stringent data selection applied in this study as discussed in section 4.2.3 may have biased the data. On a 
regional spatial scale, Baran et al. (2015) investigated the local variation of crystals roughness inferred from multi-directional 
total reflectances in relation to relative humidity. Using a rather limited dataset off the coast of Great Britain, they concluded 
that the occurrence of pristine crystal mixtures are associated with relatively humid conditions. However, these interesting con- 
clusions are based on only 12 POLDER pixels with inferred pristine particles and need to be confirmed using a larger dataset. 
In addition, biases from interference of lower liquid clouds could not be convincingly excluded. 

From the list of dominating shapes in Table | and the discussion of remote sensing studies above, it is clear that the inferred 
shapes are highly dependent on which shapes are included in the investigation. For example, Cole et al. (2014) find aggregates 
of columns to be dominating globally but do not consider Koch fractals, while the opposite is true for Chepfer et al. (2001). 
Also, Cole et al. (2013) found a rough general mixture to fit most global data and severely roughened aggregates of columns 
to be a poor fit, while Cole et al. (2014) applied a wider range of roughness to all particles and came to the exact opposite 
conclusion, i.e., a dominance of aggregates of columns and hardly any data to be more consistent with the general mixture of 
habits with any roughness applied. Furthermore, Cole et al. (2013, 2014) found the inferred shapes and roughness values to be 
dependent on assumed crystal size. However, as also discussed in section 3, this size dependency largely stems from the fact 
that most of the assumed ice crystal model geometries (i.e., aspect ratios of the crystals or their components) depend on size. 
Thus, assuming a different size essentially changes the set of models considered in this case. As discussed in section 2, the 
geometry of the components of such complex crystals, and particularly their size dependency, is very poorly constrained. A 
general aspect of the studies discussed above is that they focus on the inference of either ice crystals habit, crystal roughness or 
a combination of the two. However, as discussed in section 3, the aspect ratio of the crystals or their components greatly affects 


the scattering properties and needs to be taken account in a systematic way when inferring ice crystals shapes from remote 
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sensing data. Without the systematic inclusion of aspect ratios and roughness values in the retrieval products the results from 
such remote sensing studies are difficult to interpret as demonstrated above. 

A retrieval algorithm to infer the aspect ratios of crystals or their components in addition to particle roughness level from po- 
larized reflectances is presented by van Diedenhoven et al. (2012a) and applied to measurements of RSP (van Diedenhoven et al., 
2013) and POLDER-PARASOL (van Diedenhoven et al., 2014b). In essence, this approach uses hexagonal plates and columns 
as proxies for the components of more complex particles. Other than most other studies in Table 1, the method includes a large, 
nearly continuous range of aspect ratios and roughness values. A look up table of simulated polarized reflectances is used based 
on a database of optical properties for hexagonal plates and columns, which is calculated using the geometrics optics code de- 
veloped by Macke et al. (1996b). The aspect ratio of columns is varied between | and 50 with 26 geometrically increasing 
steps. The aspect ratios of plates are the inverse of those for columns, for a total of 51 aspect ratios. The roughness parameter, 
as defined by Macke et al. (1996b), is varied between 6 =0 and 6 =0.7 in steps of 0.05. The aspect ratio and roughness param- 
eter values that produce the simulated polarized reflectances that lead to the best fit to the measurements are considered the 
retrieved values. Furthermore, the asymmetry parameter is derived from the retrieved aspect ratio and roughness parameter per 
Fig. 3. 

The retrieval technique was evaluated by van Diedenhoven et al. (2012a) using simulated measurements based on optical 
properties of smooth, moderately roughened and severely roughened solid plates, solid and hollow columns, solid and hollow 
bullet rosettes, droxtals, aggregates of columns and aggregates of plates, as well as several mixtures of these habits (Baum et al., 
2005, 2011; Yang et al., 2015). The evaluation showed that particles with plate-like, column-like, smooth and rough compo- 
nents are generally correctly identified. For all particles the retrieved roughness parameters increase with increasing roughness 
of the particles assumed in the simulated measurements, as expected. As seen in Fig. 6, the ice crystal asymmetry parameters 
are generally retrieved to within 5%, or about 0.04 in absolute terms, largely independent of calibration errors, range and sam- 
pling density of scattering angles and random noise in the measurements. Since the asymmetry parameter is largely determined 
by the aspect ratio and roughness of the crystal components, this good agreement between retrieved and true asymmetry pa- 
rameters suggests that aspect ratio and roughness parameters of the components of these complex particles are retrieved well 
by the method, although they were not explicitly evaluated by van Diedenhoven et al. (2012a). A detailed study on simulated 
clouds consisting of bullet rosettes with arms of varying aspect ratios and roughness values showed that the algorithm retrieves 
aspect ratios with a bias of 20% on average and the roughness parameter within 0.05 when cloud optical thickness is above 5 
(van Diedenhoven et al., 2012a). Errors on retrieved aspect ratios increase with decreasing optical thickness. Interestingly, the 
method finds that the polarized reflectances of clouds consisting of hollow columns and hollow bullet rosettes most closely 
resemble those consisting of columns with high aspect ratio (a > 15). This result makes sense since the walls of the hol- 
low parts of the crystals are thin, resembling columns with high aspect ratios (see section 3 and Figs. 4 and 5). In addition, 
van Diedenhoven et al. (2012a) showed that asymmetry parameters of mixtures of smooth and rough complex particles as de- 
fined by (Baum et al., 2011) are also mostly retrieved within 5%. For mixtures of hexagonal columns and plates with varying 
aspect ratios and roughness values, van Diedenhoven et al. (2016a) showed that the average absolute errors between retrieved 


and the ensemble-average aspect ratio a<, are generally below 0.1. Furthermore, ensemble-average roughness parameters are 


30 


15 


@ plates 
asolid columns 
0.95[ Ysolid bullet rosettes 
. ¢ droxtals 


Moaggregates of columns 

Osmoll aggregates of plotes .~ A 
0.90} Alarge aggregates plates 

Vhollow columns 7 


a © hollow bullet rosettes. 
D =” 
3 0.85 
oO 
"4 

0.80 

s Ro Smooth 
0.75." ae Moderately rough 


Rough 


0.70 : 
0.70 0.75 0.80 0.85 0.90 0.95 
True g of complex ice 


Figure 6. Asymmetry parameters retrieved with the method of van Diedenhoven et al. (2012a) from simulated data based on 
optical properties of several complex crystals plotted against the true asymmetry parameters of the complex crystals. Black, 
green and red colors are results for simulated data assuming smooth (a = 0), moderately rough (a = 0.03), and severely rough 
(a = 0.5) particles, respectively. For each combination of roughness and habit, 10 different size distributions are applied. See 
van Diedenhoven et al. (2012a) for further details. The solid line shows the 1-1 line. Dotted lines indicate the targeted 5% 


accuracy limits. Figure reproduced from van Diedenhoven et al. (2012a). 


generally retrieved within 0.1. Generally, the approach tends to be somewhat biased toward retrieving column-like crystals, 
although for about 75% of the test cases the dominating geometry was correctly determined. Furthermore, only considering 
mixtures that are dominated for over 3/4 by either plates or columns yielded the correct dominating geometry in about 90% of 
the cases. 

van Diedenhoven et al. (2013) further evaluated this approach applied to measurements of the RSP instrument collected 
during the CRYSTAL-FACE campaign based in Florida in 2002. Four case studies were analyzed: two cases of thick convective 
clouds and two cases of thinner (detached) anvil cloud layers. In all cases the measurements indicate roughened ice crystals, 
consistent with previous findings. Retrieved aspect ratios in three cases were found to be close to unity, indicating that compact 
particles dominate the radiation, qualitatively consistent with CPI images where available. Retrievals for one contrasting anvil 
case indicate ice crystals consisting of plate-like components with aspect ratios around 0.3, consistent with the increased 
number of aggregates of plates seen in the CPI images obtained in this cloud layer. 

An example of crystal properties varying per cloud type and conditions is illustrated in Fig. 7, which is derived from (previ- 
ously unpublished) data of RSP and the Cloud Physics Lidar (CPL), both mounted on NASA’s high-altitude ER-2 aircraft dur- 
ing the Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys (SEAC*RS) 


campaign based out of Houston, Texas in 2013 (Toon et al., 2015). This figure also shows retrievals of ice effective radius (Eq. 
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Figure 7. Histograms of ice cloud properties retrieved from measurements of RSP and CPL obtained on 2 September, 2013, 
during the SEAC*RS campaign. From left to right respectively, the panels show the effective radii, asymmetry parameters, 
aspect ratios, roughness (or distortion) parameters 6, specular reflection index, cloud top heights, lidar penetration depths and 
depolarization ratios. Effective radii retrieved with the RSP channels at 1.59 and 2.26 zm are shown in yellow and blue, 
respectively. CPL retrieved quantities are shown in green. The top and middle panels show retrievals obtained during flight 
legs over a frontal cloud system in developing and dissipating stages, respectively. The bottom panel shows data for convective 


clouds over land. See text for more details. 


12) using the shortwave infrared measurements on RSP (Nakajima and King, 1990; van Diedenhoven et al., 2014b, 2016b) as 
well as an index quantifying the strength of observed specular reflection in the RSP data (cf. Bréon and Dubrulle, 2004). For 
completeness, the cloud top height derived by RSP (Alexandrov et al., 2012; Sinclair et al., 2017) and CPL are also given, in 
addition to the CPL lidar depolarization and penetration depth. Here, CPL penetration depth indicates the physical depth at 
which CPL signals saturate, which is related to the mean extinction at cloud top (cf. van Diedenhoven et al., 2016b). Three 
different cloud conditions were observed on September 2nd, 2013. The top row of Fig. 7 shows retrievals for a frontal cloud 
system that was sampled, indicating a rather large variability in effective radius with a peak near 35—40 jum, cloud top heights 
at around 11 km and almost exclusively rough compact crystals with aspect ratios near unity, roughness parameters around 0.5 


and asymmetry parameters near 0.75. The second row shows retrievals for the following hours, when the cloud system was 
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dissipating. Interestingly, the tops of these clouds seem to sublimate first, leaving the lower ice cloud layers visible to the RSP. 
The retrievals indicate plate-like particles with lower aspect ratios and a broader range of roughness values compared to the 
previous case. Observations of specular reflection indicates horizontally oriented ice plates were present in some regions. CPL 
depolarization ratios for these cloud layers are also lower compared to the earlier measurements, which might be caused by a 
change in ice habit, but may also be due to horizontally oriented crystals, especially in the case of depolarization values below 
0.25. RSP effective radius values for these clouds range between 15 and 45 j:m and considerable differences are seen between 
sizes retrieved with 2.26 and 1.59 jzm bands, which implies substantial vertical variations of ice sizes (van Diedenhoven et al., 
2016b). Finally, the third part of the day was devoted to sampling convection over land. The retrievals, shown in the third 
row, generally yield much smaller effective radii around 20 jzm with almost no difference between retrievals using different 
spectral bands, indicating little vertical variation. This is consistent with relatively shallow lidar penetration depths for this 
case indicating compact and opaque cloud tops. The aspect ratios and asymmetry parameters show more variation than for the 
rest of the day, although lidar depolarization at cloud top has a narrower distribution. Furthermore, roughness parameters are 
somewhat larger for these convective clouds compared to the frontal system shown in the top row. This case study demonstrates 
the complex variability that can be observed in ice shapes and other properties of ice clouds. Furthermore, it shows the benefits 
of combining different retrievals techniques and instruments crucial to obtain a more complete view of ice cloud properties. 
van Diedenhoven et al. (2014b) presented retrievals of ice crystals shape and asymmetry parameters from POLDER mea- 
surements collected off the north coast of Australia in relation with the Tropical Warm Pool - International Cloud Experiment 
(TWP-ICE) campaign in 2006 (May et al., 2008). The data are divided into periods of 4—9 days with alternating weak and 
strong convection, indicated by observed rain rates. Furthermore, the data is presented as a function of cloud top pressure and 
temperature, as ice crystal properties are generally observed to vary with temperature (e.g., Lawson et al., 2010; Noel et al., 
2004, see also section 2). Only clouds with an optical thickness larger than 5 are included. The mean results shown in Fig. 8 
indicate that mostly plate-like particle components with mean aspect ratios (w<1) around 0.6 and low asymmetry parameters 
characterize strongly convective periods, while weakly convective periods generally show particles with larger asymmetry pa- 
rameters, lower component aspect ratios, somewhat lower roughness parameters and more column-like crystal components. 
The abundance of compact plate-like crystals in the tops of convective clouds is consistent with previous observations of the 
dominance of compact and aggregated ice crystals with plate-like components observed in tropical deep convection (e.g., 
Noel et al., 2004; Connolly et al., 2005; Um and McFarquhar, 2009; Baran, 2009, see also section 2). There appears to be a 
trend towards lower aspect ratios and more column-like particles at warmer temperatures. The more column-like ice crys- 
tals with component aspect ratios further deviating from unity as indicated by the observations during the weakly convective 
periods may be consistent with a stronger contribution of particles grown in situ, which are more likely to form as crystals 
with column-like components, such as bullet rosettes, at the observed temperatures (Bailey and Hallett, 2009; Baran, 2009; 
Gallagher et al., 2012; Fridlind et al., 2016, see also Fig. 1). Comparing the results for the three strongly convective periods 
shows that microphysical parameters observed during the first two are very similar, while the third period shows somewhat 
greater roughness, fewer column-like crystals and lower asymmetry parameters. As shown by van Diedenhoven et al. (2014b), 


this later period also had substantially larger effective radii, especially at the warmer temperatures. The meteorological quantity 
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Figure 8. Mean ice crystal-component aspect ratio (a<1, top left), roughness parameter (top right) and asymmetry parameter 
(bottom right) retrieved from POLDER-PARASOL data off the north coast of Australia between 16 January and 20 February, 
2006. The percentage of retrieved column-like aspect ratios (i.e., ~@ > 1) is shown in the bottom left panel. Data within 25—hPa— 
wide cloud top pressure bins are averaged to produce profiles for five different periods (indicated by colors) with alternating 


strong and weak convective strengths. See text and van Diedenhoven et al. (2014b) for further details. 


that possibly distinguishes the third strongly convective period from the other periods with strong convection is the middle- 
to-upper tropospheric zonal wind shear, which is much weaker for the third period and may have affected crystal evolution. 
For strongly convective periods, the roughness parameter values are about 0.55 on average and decrease significantly with 
increasing cloud top temperature, while asymmetry parameters increase. These results suggest systematic variations of crystals 
shape characteristics in relation to cloud top heights and atmospheric conditions. Such relationships need to be substantiated 
with more data. As demonstrated by van Diedenhoven et al. (2014b), such observed variations have significant impacts on 
the radiative properties of convective clouds, which need to be better understood to improve their representation in climate 
predictions. 

In order to compare convective clouds observed in different atmospheric regimes, Fig. 9 shows histograms of ice crystal 


asymmetry parameters, aspect ratios and roughness parameters retrieved over convective clouds using 1) POLDER-PARASOL 
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Figure 9. Histograms of the retrieved asymmetry parameters (left), aspect ratios (middle) and roughness parameters (right) for 


convective clouds observed in relation to the TWP-ICE (red), TC* (green) and SEAC*RS (yellow) campaigns. 


data collected in relation with the TWP-ICE campaign (May et al., 2008; van Diedenhoven et al., 2014b); 2) POLDER-PARASOL 


data collected in relation to the TC* campaign (Toon et al., 2010) in 2007 near Costa Rica (previously unpublished); and 3) RSP 
data collected during the SEAC*RS campaign (Toon et al., 2015) based out of Houston, Texas in 2002 (cf. van Diedenhoven et al., 
2016b). Only ice clouds with optical thicknesses larger than 5 are included in all data sets. Statistics of asymmetry parame- 
ters, aspect ratios and roughness parameters derived from these different datasets are very consistent. The data indicate crys- 
tals with mostly plate-like components (a@ < 1), which is consistent with in situ measurements in convective clouds (e.g., 
Um and McFarquhar, 2009). Also, aspect ratios close to unity are mostly found, indicating the dominance of compact parti- 
cles. Roughness values are generally large with maxima greater at 0.5, which is largely consistent with the roughness statistics 
found in tropical regions by Cole et al. (2014). Asymmetry parameters peak at around 0.76, but the distribution shows a sub- 
stantial tail toward larger values. 

In summary, results from retrieval approaches that focus on retrieving specific ice crystal habits are generally inconclusive 
and highly dependent on the shapes included in investigations. One general conclusion that can be derived from the various 
studies is that particle roughness is prevalent. However, roughness levels are found to vary with location, cloud top temperature, 
and atmospheric conditions. Systematic retrievals of crystal component aspect ratios show convective cloud regimes generally 
have crystals with aspect ratios close to unity and mostly plate-like particles, but the particle properties depend on multiple 
factors, e.g., cloud type, cloud height, convective strength and possibly other dynamical quantities as wind shear, as well as 
humidity. More global and local studies are needed to untangle such relationships between ice crystals shape and cloud type, 


cloud height and atmospheric conditions. 
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5 Prospective 


Although many studies on remote sensing of ice crystals shapes using lidar and multi-directional reflectance data have been 
performed over the past few decades, the discussions in sections 4.1 and 4.2 show that it remains difficult to extract systematic 
conclusions from these studies. Specifically, robust quantitative interpretations of lidar depolarization measurements are still 
not available. As discussed in section 4.1, more research on optical properties calculations will likely improve the prospective 
of gaining quantitative information about ice crystals shapes from lidar measurements. Furthermore, measurements of a multi- 
Static lidar, as proposed by Mishchenko et al. (2016), measuring the backscattered signal at additional angles and thus probing 
the depolarization properties at two or more angles could yield increased information content for the retrieval of ice crystal 
shape from lidars (cf. Smith et al., 2016). Combining colocated lidars at multiple wavelengths is not expected to increase 
potential for ice shape retrievals since the ice scattering properties are largely wavelength independent. Also high spectral 
resolution capability, such as included in ATLID and the airborne HSRL (Burton et al., 2015), is not expected to increase 
information content on ice crystal shapes, although it will provide valuable measurements of ice cloud extinction. 

Currently, the POLDER instruments have been the only multi-directional polarimeters deployed in space. As polarimetry 
has great potential for cloud retrievals, as well as for the inference of aerosol properties, polarimeters are considered for many 
future satellite missions. For instance, the Multi-viewing, Multi-Channel Multi-Polarization Imaging instrument (3MD) is a 
follow-up version of POLDER and is planned to be included on the European MetOp series (Marbach et al., 2013). NASA’s 
upcoming Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission will likely include a multi-directional polarimeter. 
Furthermore, the Hyper-Angular Rainbow Polarimeter (HARP) instrument is a cubesat mission that is slated be be launched 
soon. Also the Multi-Angle Imager for Aerosols (MAJA) instrument is selected by NASA for further development and space 
deployment. It is expected that all these multi-directional polarimeters have potential for the retrieval of ice crystal shapes that 
is similar to or better than POLDER’s. Interestingly, the 3MI instruments will in include a channel at around 1370 nm that 
is located on a strong water vapor absorption band. This limits the band’s sensitivity to the surface and to liquid clouds in 
the lower atmosphere and increases the sensitivity to thin cirrus (cf. Gao et al., 1993). The RSP band at 1880 nm has similar 
capabilities and was used by Ottaviani et al. (2012) to infer the ice shape of a thin cirrus layer over the Deepwater Horizon 
oil spill site in 2010. Furthermore, the SPEX airborne polarimeter (Rietjens et al., 2015) yields multi-directional polarization 
measurements in the oxygen A-band that similarly shields the lower atmosphere and surface. Further statistical evaluation of 
such measurements will provide valuable statistics of thinner cirrus clouds, which are largely excluded from the current remote 
sensing results. 

Ice crystal shape retrieval approaches have been largely limited to finding best fits to the measurements within look up 
tables of simulated measurements. It is advisable that more systematic inversion techniques are employed since these allow 
faster data processing, non-discrete solutions, and, more importantly, better error estimations (Rodgers, 2000). For example, 
Hioki et al. (2016) recently developed an algorithm based on empirical orthogonal function analysis of polarized reflectances to 
infer ice crystal roughness levels from such measurements. Approaches that retrieve quantifiable parameters such as ice crystal 


component aspect ratios and roughness parameters are especially suitable for the implementation in an algorithm employing 
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an inversion technique. The results of van Diedenhoven et al. (2016b) suggest that such an approach could infer the relative 
contribution of plate-like and columnar components to the ice crystal distributions, in addition to the mean aspect ratio (a<1) 
and roughness parameter. 

To date, few systematic investigations on the variation of crystal shape, aspect ratios and crystals roughness and their re- 
lation with cloud type, height and atmospheric conditions have been performed (e.g., Baran et al., 2016; Cole et al., 2014; 
van Diedenhoven et al., 2014b). More global and local studies are needed to untangle such relationships between ice crystals 
shape and cloud type, cloud height and atmospheric conditions. Such relationships will provide observational constraints for 
improved parameterizations of ice cloud optical properties (e.g., Baran et al., 2016) and for microphysical packages for cloud 


simulations, especially those explicitly prediction ice crystal shapes (e.g., Hashino et al., 2007, 2011; Harrington et al., 2013). 


6 Conclusions 


Improved constraints on the natural variation of ice crystal shapes is important since the shape of ice greatly affects their 
radiative and microphysical properties. Theoretical and laboratory studies show that ice crystal shapes largely depend on 
temperature and humidity of the environment in which they grew. However, in situ measurements in real ice clouds generally 
show complex mixtures of shapes and large contributions of irregular, complex, aggregated crystals. Also, the aspect ratios 
of components of these crystals are found to vary substantially. In addition, high magnification imaging of ice crystals show 
roughness structures of various forms and levels on the ice surfaces of growing and sublimating ice crystals. All of these macro- 
and microscale ice shape characteristics substantially affect the radiative properties of ice clouds and better constraints on how 
these shape characteristics vary with cloud type, temperature, humidity, locations, availability of aerosols, etcetera, are crucial 
to improve the representation of ice clouds in climate projections. This chapter reviews the current state of remote sensing of 
ice crystal macro- and microscale structure. 

The radiative properties of clouds are determined by the number of ice crystals, their extinction cross sections, single scatter- 
ing albedos and scattering phase matrices. Of these optical properties, the scattering phase matrix is especially relevant for the 
remote sensing of ice shapes, since it is substantially dependent to ice crystal shape but relatively independent to the size of ice 
crystals at non-absorbing wavelengths. The angular features in the scattering phase function and the linear polarization phase 
function that depend on ice crystals shape are used by remote sensing studies. Also the depolarization of backscattered polar- 
ized lidar signals is used to obtain information on ice crystal shape. Reviewing the dependencies of scattering phase matrices 
on ice shape reveals that the phase matrices are mostly determined by the aspect ratios of components of complex crystals as 
well as by the microscale structures such as crystal roughness, while the macroscale shape (i.e., habit) is of lesser importance. 
Furthermore, while ice particle macroscale shape has a seemingly endless variability and is not a quantifiable parameter, crystal 
component aspect ratio and roughness level are quantifiable and can be systematically related to variations in the phase matrix. 
For example, phase function asymmetry parameters increase as aspect ratios deviate from unity and decrease as roughness 
levels increase. It is therefore advised that remote sensing studies focus on retrieving information about crystal component 


aspect ratios and microscale structure, rather than on inferring the occurrence of specific ice habits. 
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A review of the literature on lidar depolarization measured in ice clouds shows that depolarization is generally found to 
increase with increasing cloud height and also varies with latitude. This variation is generally linked to the variation of ice 
crystal shape. However, the interpretation of the depolarization remains largely qualitative and inconclusive. For simple, smooth 
hexagonal prisms, lidar depolarization is shown to vary with crystal aspect ratio, but studies on the effects of crystal macroscale 
complexity, microscale roughness and hollowness are lacking, contradictory or inconclusive. More research on the relation 
between lidar depolarization and ice crystal shape is advised in order to move toward more quantitative inferences of ice 
crystal shape properties from lidar measurements. 

Numerous studies evaluated the angular variation of total and/or polarized reflectances of ice clouds in order to infer in- 
formation about ice crystal shape from them. A general conclusion is that ice crystal surface roughness or crystal distortion 
is prevalent. However, the conclusions about the dominating ice shapes are often contradictory. Furthermore, the data sug- 
gests that ice crystals shape and roughness is highly variable. Perusing the various studies, it is clear that the inferred shapes 
are strongly dependent on which shapes are included in the investigation and that such selections of shapes can be inter- 
preted as rather arbitrary. Moreover, the specific geometries to define these shapes, most importantly the aspect ratios of their 
components, are very unconstrained. Since the inferred crystal shapes often do not agree with what is expected from in situ 
measurements and theory, the retrieved shapes are often interpreted as merely "radiative equivalent effective shapes" with no 
direct relation to ice crystal shapes in the observed clouds. However, focusing remote sensing applications on retrieving crystal 
component aspect ratio and particle roughness, rather than ice habit, yields useful physical information about these shape char- 
acteristics. Retrieval approaches that focus on inferring aspect ratios and/or roughness reveal that ice shapes depend on cloud 
height, latitude, cloud type, convective strength and possibly on humidity and dynamical quantities as wind shear. Statistics of 
ice crystal component aspect ratios, roughness parameters and asymmetry parameters of convective clouds in several different 
regimes are very comparable, suggesting generalized conclusions could be derived from such measurements. The data for 
these convective clouds indicate crystals with mostly compact plate-like components with high roughness levels are prevalent. 
Asymmetry parameters peak at around 0.76, but the distribution is substantially wide with a tail toward large values. 

As is generally the case with all remote sensing products, the results in the current literature represents a subset of all 
ice clouds. Especially studies using multi-directional measurements are often biased to relatively optically thick ice clouds 
and thus excludes most cirrus. Improved data selection (e.g., Wang et al., 2014) and analysis of measurements obtained at 
wavelengths with substantial atmospheric gas absorption (e.g., Ottaviani et al., 2012) could yield more information about thin 
cirrus. Furthermore, the multi-directional total and polarized reflectances as well as lidar depolarization only yield information 
of the top 1-3 optical depths of clouds, while ice crystal shape is likely to vary vertically in clouds. This limitation also applies 
to retrievals of other cloud particle properties such as phase and size (Platnick, 2000; van Diedenhoven et al., 2016b). Hence, 
it is important to augment such retrievals with, e.g., in situ studies of the vertical variation of cloud particle properties. For 
convective clouds, however, statistics of cloud top trends of ice crystal properties can generally be used as surrogates for 
trends with height within convective cloud tops (referred to as time-space interchangeability, Lensky and Rosenfeld, 2006; 
van Diedenhoven et al., 2016b). Furthermore, cloud top properties are very relevant as the top of clouds can be considered as 


radiatively the most relevant part. 
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Lidar depolarization and multi-angular measurements, in addition to in situ observations, consistently reveal that ice crystal 
shape varies considerable with cloud top height, cloud type, location and atmospheric state. Ignoring these variations leads to 
biases in retrievals of ice effective radius and optical thickness that are sensitive to shape as well as to biases in simulations 
of ice cloud properties and their radiative effects. Research has begun to untangle such relationships between ice crystals 
shape and cloud type, cloud height and atmospheric conditions, but more studies are needed to reach systematic conclusions. 
Given the high potential of polarized lidar and multi-angle polarimeters for cloud and aerosol remote sensing, such instruments 
are planned or considered for many future satellite missions providing many future opportunities to further study the global 


variation of ice crystal shapes. 
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